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Abstract 


\  t 

The  forecasting  of  the  performance  of  electric  vehicles  has 
revealed  many  discrepancies  when  compared  to  actual  test  data.  An 
electric  vehicle  computer  simulation  program  was  developed  to  ameliorate 
this  deficiency.  The  approach  was  to  establish  a  very  comprehensive  and 
flexible  vehicle  model  and  simulate  its  operation  on  a  realistic  driving 
cycle.  The  driving  cycle  selected  was  the  Federal  Urban  Driving  Sequence. 
A  thorough  vehicle  model  was  established  that  incorporates  aerodynamic 
drag,  rolling  resistance,  both  rotational  and  translational  inertial 
effects,  and  component  by  component  dynamic  power  train  efficiencies. 
Battery  discharge  performance  is  tracked  by  a  fractional-utilization 
algorithm  with  corrections  for  short-term  discharge  effects.  The  simu¬ 
lation  compares  required  power  obtained  from  the  driving  cycle  speed 
schedule  and  vehicle  model  characteristics  with  the  available  power  at  the 
motor  for  each  time  increment  of  the  driving  cycle  to  compute  battery 
fraction  used  and  deviation  from  the  speed  schedule  when  available  power 
is  insufficient.  These  results  of  the  simulation  can  be  used  to  evaluate 
an  existing  vehicle's  performance  or,  if  desired,  to  develop  vehicle 
parameters  to  obtain  a  specified  performance  level.  An  application  of  the 
program  to  develop  a  suburban  passenger  vehicle  is  included  to  demon¬ 
strate  the  simulation's  utility.  A  test  bed  vehicle  was  constructed  and 
tested  to  verify  the  simulation.  Additional  aspects  such  as  microprocessor 
based  controllers  including  implementation  of  an  optimal  control  law  were 
investigated  to  gain  insite  into  the  efficiency  and  performance  trade-offs 
of  such  a  system. 
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Va 

- 
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vc 

= 
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Vm 

= 
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VM 
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Vehicle  mass  (for  test  conditions)  (kg) 
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Energy  (watt  *  sec) 
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ELECTRIC  VEHICLE  MODELING 
AND  SIMULATION 


I .  Introduction 

This  report  describes  recent  efforts  of  the  Air  Force  Institute  of 
Technology  in  the  area  of  Electric  Vehicle  (EV)  research  aimed  at  developing 
models  to  aid  government  agencies  in  petroleum  fuel  conservation  for  ground 
transportation.  Interest  in  petroleum  fuel  conversation  and  alternative 
fuel  usage  has  greatly  increased  in  the  past  decade.  The  USAF  and  other 
government  agencies  are  becoming  especially  interested  in  reducing  foreign 
oil  dependence.  While  NASA  and  USAF-sponsored  projects  have  achieved 
significant  technological  advances  toward  reducing  aircraft  fuel  usage, 
only  recently  has  complimentary  interest  been  shown  toward  economic 
alternatives  to  the  gas-powered  ground  vehicle. 

Background 


In  the  mid-1970's,  replacement  of  various  small  conventional 
vehicles  with  electric-driven  equivalents  became  an  item  of  interest  with 
Air  Force  Logistics  Command  (AFLC/USAF).  The  maintenance  engineering 
branch  of  AFLC  conducted  several  projects  aimed  at  reducing  operating 
costs  of  light  intra-base  utility  vehicles  by  substitution  of  standard 
size  vehicles  with  low-powered  gasoline  or  electric  vehicles  (Ref  33). 
General  conclusions  were  favorable  toward  the  small  electric  vehicles 
even  though  their  operating  costs  per  mile  during  one  14-month  test 
period  were  actually  somewhat  higher  than  the  gas-powered  alternatives 
(Ref  6).  Further  evidence  of  USAF  electric  vehicle  interest  came  in  May 
1978,  when  an  Air  Command  and  Staff  College  research  study  investigated 
the  practical  and  economical  feasibility  of  using  these  vehicles  for 
several  on-base  missions  (Ref  34).  This  study  concluded  that  electric 
vehicle  operating  costs  were  nearly  equal  to  standard  vehicles,  but  that 
the  higher  initial  investment  and  limited  performance  of  electric  vehicles, 
due  partly  to  their  rudimentary  state-of-the-art  development,  made  them 
undesirable  at  that  time.  (Note  that  the  operating  expenses  for  the 
above  study  were  based  on  energy  costs  of  4c  per  kilowatt  hour  of 
electricity  and  40c  per  gallon  of  gasoline).  A  further  recommendation  of 
the  report,  however,  was  that  the  USAF  become  actively  involved  in  the 
electric  vehicle  tests  sponsored  by  the  Department  of  Energy  (DOE). 

This  DOE  program  was  established  by  Public  Law  94-413,  "Electric 
and  Hybrid  Vehicle  Research,  Development,  and  Demonstration  Act  of  1976." 
Among  other  things,  it  authorized  the  Energy  Research  and  Development 
Agency  (ERDA,  now  absorbed  into  DOE)  to  test  7,500  electric  and  hybrid 
vehicles  within  five  years.  Public  Law  95-238,  signed  in  February  1978, 
amended  this  Act,  including  an  extension  of  the  demonstration  phase 
through  1986.  In  July  1980,  by  DOE  interagency  Agreement  DE-AIOI- 
80CS50208  (Ref  8),  the  USAF  became  formally  involved  with  the  demonstration 
program.  This  agreement  committed  the  USAF  to  procure  and  test  15  electric 
vehicles  to  be  distributed  to  three  different  locations  for  a  total  period 
of  four  years. 


An  early  product  of  Public  Law  94-413  tasking  was  a  very  thorough 
NASA  report  titled,  "State-of-the-Ar t  Assessment  of  Electric  and  Hybrid 
Vehicles"  published  in  September  1977.  The  purpose  of  this  report  was  to 
characterize  the  current  state  of  electric  and  hybrid  vehicle  development 
through  controlled  tests  of  commercially  available  and  experimental 
vehicles,  through  manufacturer-supplied  data,  and  through  user  reports. 
Data  was  thus  obtained  from  nearly  200  different  electric  and  hybrid 
vehicles.  Of  particular  significance,  the  report  discovered  extreme 
variances  between  similar  vehicles  in  their  overall  energy  efficiency 
(.015  to  .8  kwh/km/1000  kg)  and  in  their  user-experienced  performance 
versus  manufacturers'  claims  or  controlled  laboratory  test  results. 

Objectives 


In  lieu  of  the  increased  interest  among  government  agencies 
(including  very  recently  the  USAF)  as  well  as  the  observed  lack  of 
sophistication  among  many  electric  vehicle  designs,  the  prime  objective 
of  this  effort  was  to  provide  an  accurate  method  to  simulate  electric 
vehicle  systems.  The  goal  was  to  provide  a  computer  simulation  useful 
both  for  new  system  design  as  well  as  evaluation  of  existing  competitive 
designs.  To  offer  these  capabilities,  the  model  had  to  be  flexible  enough 
to  accommodate  a  wide  variety  of  design  options. 

Preliminary  results  of  the  simulation  were  used  to  select  a  viable 
electric  vehicle  system  to  compete  economically  with  conventional  USAF 
passenger  cars.  This  system  was  then  incorporated  into  a  low-cost 
vehicle  which  was  designed  and  built  for  data-gathering  and  simulation 
validation  as  well  as  a  flexible  test-bed  to  investigate  other  EV  options 
such  as  controller  performance. 

Approach 

Based  upon  the  lack  of  correlation  between  actual  performance  and 
early  simulation  prediction  attempts  (Ref  32),  an  early  priority  of  this 
effort  was  to  carefully  select  an  operating  sequence  to  most  realistically 
portray  on-the-road  vehicle  use.  Once  determined,  road  load  power  require¬ 
ments  could  be  easily  calculated.  The  dynamic  efficiency  of  each  power 
train  component  was  then  individually  modeled  and  assembled  together  into 
an  interrelated  complex  vehicle  model.  Net  power  requirements  to  the  motor 
and/or  engine  could  then  be  obtained.  A  multi-variable  battery  model  was 
combined  with  controller  characteristics  to  determine  available  power  to 
the  motor  under  any  condition.  A  comparison  between  this  available  power 
and  the  required  power  was  used  to  predict  the  sufficiency  of  the  power 
supply  and/or  the  performance  deviation  from  the  prescribed  driving  cycle 
(speed  schedule).  Further  features  of  the  battery  model  enabled  the 
simulation  to  track  net  energy  usage  and  the  extent  of  nominal  discharge. 
Final  outputs  of  the  simulation  were  then  the  average  energy  cost  per 
kilometer  as  well  as  the  amount  of  battery  consumed  per  cycle,  which 
together  with  estimated  procurement  and  maintenance  costs,  can  be  used 
as  a  measure  of  merit  for  a  selected  sv  tern. 


The  basic  approach  toward  achieving  the  objective  of  test 
vehicle  design  and  construction  **as  to  use  an  existing  chassis  combined 


with  readily-available  electric  power  components.  This  approach,  while 
not  providing  a  completely  system-optimized  vehicle,  allowed  quick 
construction  of  a  low-cost  test-bed  for  further  component  optimization 
investigation. 

Performance  Requirements 

The  flexibility  of  the  vehicle  model  and  simulation  developed  later 
in  this  study  allows  their  use  for  a  variety  of  vehicle  types  and  missions. 
However,  the  prime  vehicle  envisioned  throughout  the  study  and  which  formed 
the  basis  for  the  test  vehicle  was  a  compact  four-passenger  sedan.  This 
type  of  general  purpose  vehicle  is  used  by  the  USAF  for  both  on-  and  off- 
base  missions.  Table  I  lists  the  performance  goals  that  were  established 
for  this  type  of  vehicle  based  on  the  expressed  desires  of  USAF  vehicle 
managers  (Ref  34)  and  analysis  of  some  reasonable  off-base  mission  require¬ 
ments.  A  vehicle  designed  to  minimally  meet  the  Table  I  performance  would 
be  primarily  an  urban  or  on-base  vehicle,  but  would  be  capable  of  certain 
high  speed,  extended  range  trips.  The  90  km/hr  top  speed  with  80  km/hr 
cruise  allows  limited  use  of  the  vehicle  on  typical  suburban  routes,  while 
the  0  to  50  km/hr  acceleration  satisfies  typical  urban  driving  demands 
(Ref  32). 

TABLE  I 


Baseline  Performance 


Payload 

4  Passenger  (270  kg) 

Top  Speed 

90  km/hr 

Acceleration 

i  0  to  50  km/hr 

15  sec 

Range 

Urban 

80  km 

80  km/hr  cruise 

60  km 

For  a  test  vehicle  with  an  estimated  gross  mass  of  1200  kg  and 
frontal  area  of  1.9  m2,  published  data  (Refs  17,  23)  indicated  that  the 
90  km/hr  speed  translated  to  approximately  12  kw  of  required  motor 
power,  whereas  the  0  to  50  acceleration  required  approximately  22  kw. 
The  apparent  critical  parameter  in  electric  vehicles,  however,  is  total 
energy  storage  capacity  rather  than  maximum  power  available  (Ref  15). 
Multiplying  the  two-hour  high  speed  cruise  requirement  by  an  estimated 
10.5  kw  power  requirements  for  80  km/hr  resulted  in  a  net  required 
energy  storage  of  21  kwh.  These  power  and  energy  requisites  based  on 
the  Table  I  performance  then  became  important  criteria  for  selection  of 
the  electric  drive  system  for  the  test  vehicle. 

With  the  vehicle  performance  requirements  defined,  the  modeling 
of  appropriate  classes  of  subsystems  could  now  proceed. 


II.  Modeling 


The  basis  for  any  simulation  is  established  by  the  models  it 
incorporates.  The  EV  simulation  is  composed  of  models  of  vehicle  retarding 
forces,  vehicle  body  and  subsystems  such  as  drive  train,  controller, 
batteries  and  hybrid  engine. 

Derivation  of  Road  Loads 

For  travel  on  level  roads,  the  combined  vehicle  retarding  forces 
can  be  divided  into  the  rolling  resistance,  air  resistance,  and  acceleration 
forces.  Road  grade  and  wind  were  not  explicitly  modeled  in  this  effort 
because  standard  driving  cycle  tests  do  not  include  them.  In  Chapter  V, 
both  grades  and  wind  gusts  were  included  as  disturbances  when  investigating 
controller  robustness. 

Rolling  Resistance.  The  force  required  to  move  any  vehicle  at 
constant  speed  without  considering  aerodynamic  or  powertrain  drag  is 
developed  from  the  work  performed  in  flexing  the  tire  casing,  the  work  in 
compressing  the  road  surface  under  the  tire,  and  the  frictional  work  from 
tire  slip.  As  such,  this  retardant  force  is  logically  a  function  of  tire 
design,  inflation  pressure,  vehicle  weight,  rolling  surface,  and  vehicle 
speed.  The  dimensionless  parameter  that  combines  these  variables  is  the 
coefficient  of  rolling  resistance  CR.  The  rolling  resistance  FR  is  then 
simply: 


fR  “  CR  '  VM  •  g  (  1  ) 

where  VM  is  the  vehicle  test  mass.  Data  shown  on  Fig  1  represents  the  CR 
of  radial  ply  tires  rolling  on  a  hard,  smooth  surface  (Ref  32).  This 
information  was  used  to  obtain  the  following  continuous  equation  for  CR 
versus  speed  v: 

CR  =  CR  +  ( 2 *09  '  10_7)v2*8  (  2  ) 

Cp  represents  the  limiting  value  of  CR  at  zero  speed  (in  the  case  of  Fig  1, 
CR  would  equal  .0105).  The  speed  v  in  Eq  (  2  )  is  expressed  in  m/sec. 
Combining  Eqs  (  1  )  and  (  2  )  then  results  in  the  final  equation  for 
rolling  resistance: 


Fr  =  {C&  +  (2.09  •  10"7)v2'8}VM  •  g  (  3  ) 

which  results  in  newtons  force  with  VM  expressed  in  kilograms.  Strictly 
speaking,  Eq  (  3  )  is  correct  only  for  straight-ahead  driving;  during  ar.y 
turn  the  tires  must  also  generate  a  side  force  through  a  slip  angle. 

This  side  force  then  combines  vectorially  with  FR  to  produce  the  net 
resistance.  However,  most  urban  turns  require  slip  angles  of  less  than 
five  degrees  which  produces  a  negligible  increase  in  resistance.  While 
CR  is  left  variable  to  accommodate  different  tire  types,  it  also  can  be 
modified  to  include  grade  resistance.  Since  grade  resistance  is  simply 
the  downhill  component  of  the  vehicle  weight,  i.e.  the  sine  of  the  slope 
times  the  vehicle  weight,  the  sine  of  the  slope  may  be  added  to  CR  to 
arrive  at  the  combined  rolling  resistance. 


Aerodynamic  Drag.  At  constant  speeds  above  60  km/hr,  air 
resistance  becomes  the  dominant  retarding  force  for  most  vehicles  (Ref  28). 
This  source  of  drag  comes  from  both  viscous  skin  friction  as  well  as  form 
or  pressure  drag.  The  form  drag  caused  by  the  low  pressure  zone  behind 
the  vehicle  is  dominant  factor  contributing  approximately  90  percent  of 
the  total  airdrag  (Ref  42).  In  linearized  form,  the  total  air  resistance 
force  Fd  is  given  by: 


F  =  (£■  2 

*D  *2)  v 


(  4  ) 


where  p  is  the  ambient  air  density,  A  is  the  projected  vehicle  area  in  the 
direction  of  motion,  and  is  a  dimensionless  parameter  accounting  for  the 
shape  of  the  vehicle.  The  air  density  p  is,  of  course,  a  function  of 
ambient  pressure  and  temperature.  However,  normal  atmospheric  pressure 
variance  causes  only  minimal  density  changes.  Therefore,  replacing  p 
in  Eq  (  4  )  with  its  value  at  standard  pressure,  while  maintaining  its 
temperature  dependence,  results  in: 


F 


D 


,  176.4 
''273TrTT 


(  5  ) 


which  provides  newtons  force  with  temperature  T  in  Celsius. 


Acceleration  Forces.  In  realistic  urban  driving  cycles,  the 
inertial  resistance  to  acceleration  will,  at  times,  be  much  larger  than 
either  the  rolling  or  air  resistance.  Besides  the  translational  effects 
on  the  vehicle  mass: 


F 


AC 


VM(£) 


(  6  ) 


the  angular  accelerations  applied  to  the  rotating  masses  also  produce 
significant  inertial  reactions.  The  acceleration  torque  due  to  rotary 
inertia  is: 


M 


AC 


= 


(  7  ) 


where  J  is  the  moment  of  inertia  for  each  rotating  mass  and  —  is  the 
respective  angular  acceleration.  Chapter  IV  goes  into  more  depth  on  the 
modeling  of  this  rotary  inertia  term  and  how  its  effect  is  combined  with 
Eq  (  6  ). 


Calculation  of  Road  Load  Power.  For  use  in  the  simulation,  the 
calculated  resistive  forces  (Eqs  3,  5,  and  6)  need  to  be  converted  into  a 
required  power  as  a  function  of  time.  An  incremental  approach  was  used 
which,  for  simplicity,  assumes  constant  acceleration  for  each  one-second 
time  interval.  The  average  speed  can  then  be  used  in  Eqs  (  3  )  and  (  5  ) 
while  the  acceleration  term  in  Eq  (  6  )  becomes  simply  the  incremental 
speed  difference  divided  by  the  one-second  time  interval.  Thus,  the 
combined  road  forces  for  any  increment  are  reduced  to  an  average  force 
over  that  increment.  Required  power  can  then  be  obtained  by  multiplying 
this  combined  force  times  the  average  velocity.  This  mid-interval  power 
is  not  assumed  constant  throughout  the  increment,  thus  the  input  energy 
W  is  calculated  by: 


6 


W  =  /Pf’t 


(  8  ) 


where  the  integration  of  power  P  is  taken  over  each  time  increment. 
Substituting  P  =  F  •  v  where  F  is  the  combined  resistive  forces  and 

i  /  “V .  ■  _  *■*  _  r  \  ■  i  j 


v  =  t(~)  +  vo  into  Eq  (  10  )  yields: 


Av 


"  ■  +  Vdt 


(  9  ) 


F  is  assumed  constant  for  any  time  increment,  thus  integrating 


V  “V 

W  =  F^2A~t~'0)  At2  +  VoAt} 


(  10  ) 


or 


vl  +v 

W  =  F{— -2— °}t  =  F  •  v  •  At 


(  11  ) 


and  recognizing  the  P  =  F  •  v,  where  v  is  the  averge  speed  for  any  increment 
results  finally  in 


W  =  P  •  At 


(  12  ) 


This  work  W  then  represents  the  incremental  energy  required  at  the  rear 
wheels  to  meet  the  road  load  demands.  Modeling  of  the  vehicle  resistance 
parameters  and  power  train  losses  is  now  required  to  determine  the  necessary 
battery  output  to  overcome  this  road  load. 


Vehicle  Mass  and  Frontal  Area 


The  vehicle  mass  and  frontal  area  are  the  primary  parameters  used 
in  determining  road  loads.  As  shown  in  Eqs  (  3  )  and  (  6  ),  the  vehicle 
mass  effects  both  the  rolling  resistance  and  the  acceleration  loads. 

The  gross  vehicle  mass  may  be  considered  to  be  represented  by 

VM-CM+FM+PM+BM+MM  (  13  ) 


where 


GVM  =  gross  vehicle  mass 
CM  =  chassis  mass 
FM  =  fixed  mass 
PM  =  payload  mass 
BM  =  battery  mass 

MM  =  motor  mass  (to  include  controls  and  hybrid  engine  mass 
if  applicable) 

The  chassis  mass  includes  all  the  basic  structural  materials  in 
the  vehicle  which  are  designed  to  support  the  entire  gross  vehicle  weight. 
Items  such  as  the  frame,  suspension,  wheels,  and  tires  are  included  in  this 
category.  This  mass  is,  therefore,  a  function  of  the  gross  vehicle  mass 
and  for  passenger  vehicles  is  typically  (Ref  2)  related  as 

CM  =  ,23(GVM)  (  14  ) 
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However,  in  this  model,  the  chasis  mass  includes  an  adjustment  for  a 
fixed-ratio  transmission  and,  hence,  was  set  to 


CM  =  . 26(GVM)  (  15  ) 

The  fixed  mass  figure  represents  all  non-structural  vehicle 
components  which  are  related  to  the  selected  vehicle  type  and  payload 
rather  than  directly  to  the  power  system  or  gross  vehicle  mass.  Passenger 
seats,  instruments,  body  panels  and  accessories  are  included  in  this 
category. 

While  maximum  payload  mass  must  be  used  in  determining  the  gross 
vehicle  mass  and,  hence,  the  chassis  mass,  the  performance  predicted 
through  this  simulation  is  based  on  a  more  normal  one-half  maximum 
payload.  Combining  this  factor  with  Eq  (  13  )  results  in  a  vehicle  test 
mass  VM  as 


VM  =  GVM  -  PM/2  (  16  ) 

where,  from  Eqs  (  13  )  and  (  15  ) 

GVM  =  ,26(GVM)  +  FM  +  PM  +  BM  +  MM  (  17  ) 


or  finally 


VM  =  (FM  +  BM  +  MM)/. 74  +  .85(PM)  (  18  ) 

While  the  above  test  mass  satisfies  the  requirements  for  obtaining 
rolling  resistance,  and  additional  factor  of  effective  rotary  mass  must  be 
taken  into  account  to  predict  acceleration  loads.  Rotating  mechanisms  can 
conveniently  be  divided  into  two  categories:  those  which  rotate  at  speeds 
fixed  to  the  vehicle  speed  (such  as  tires,  wheels,  and  axles)  and  those 
which  rotate  at  the  engine  or  motor  speed  (such  as  clutch,  flywheel,  and 
crankshaft  or  armature).  While  it  is  possible  to  analytically  obtain 
values  for  both  categories  by  summing  individual  component  angular 
momenta  (Ref  26),  the  most  reliable  results  are  obtained  experimentally. 

The  net  effect  of  both  categories  can  be  resolved  into  an  equivalent 
mass  located  at  the  rim  of  the  wheels.  Since  this  equivalent  mass  is 
now  traveling  at  the  vehicle  speed,  it  can  be  considered  to  undergo  the  same 
acceleration  as  the  rest  of  the  vehicle  mass  does.  The  rotating  power 
train  component  mass  is  proportional  to  the  vehicle  mass  and  can,  therefore, 
be  considered  to  be  a  multiplicative  mass  correction  factor.  This  factor 
relates  shaft  speed  to  wheel  rim  speed  via  the  rolling  radius  and  overall 
gear  reduction  ratio.  Experimentally  derived  data  (Ref  28)  for  four-wheel 
passenger  vehicles  suggests  an  equivalent  inertia  effect  of  about  two 
percent  of  the  vehicle's  mass  for  the  first  category  of  wheel-rotation 
speed  components.  The  second  category  can  be  shown  to  contribute  an 
equivalent  mass,  varying  quadratically  with  the  shaft  speed/vehicle  speed 
ratio.  The  equivalent  vehicle  mass  then  to  be  considered  during  any 
acceleration  period  is  the  product  of  the  vehicle  test  mass  VM  and  the 
mass  correction  factor  MCF  where 


MCF  =  1.02  +  ,Q00294(A— — R-)2 


(  19  ) 


where  AR  and  GR  are  axle  ratio  and  transmission  ratio  respectively,  and 
r  is  the  tire  rolling  radius  in  meters.  While  the  data  (Ref  28)  to 
derive  Eq  (  19  )  above  was  obtained  from  conventional  gas-engined  vehicles, 
the  largest  contributor  to  the  second  half  of  this  equation  is  the  clutch- 
flywheel  assembly  which  is  generally  left  intact  for  multi-ratio  electric 
vehicles . 


The  frontal  area  of  any  optimally-designed  vehicle  will  also,  in 
general,  be  a  function  of  the  vehicle  mass.  While  geometric  up-scaling 
of  any  volume  (i.e.  mass)  suggests  a  corresponding  two-thirds  power 
increase  in  any  exposed  areas,  in  actuality,  more  massive  vehicles  tend 
toward  a  higher  proportion  of  non-payload  to  payload  volume  and  are,  there¬ 
fore  generally  denser.  This  tendency,  coupled  with  practical  passenger 
space/vehicle  width  considerations  which  result  in  larger  fineness  ratios 
with  more  massive  vehicles,  causes  actual  frontal  area  ratios  to  vary  at 
a  somewhat  lower  than  two-thirds  power  of  the  mass  ratios.  An  inspection 
of  14  electric  vehicles  detailed  in  Ref  32  indicates  that  approximately 


,FA\  ,VM'  1/2 

'fa  vvm 


(  20  ) 


or 


»*•  ■ 


(  21  ) 


where  FA'  represents  the  changed  frontal  area  caused  by  a  corresponding 
increase  or  decrease  in  vehicle  mass  VH'. 


Axle  and  Transmission  Characteristics 


Energy  losses  in  the  conventional  automotive  rear  axle  occur  from 
relatively  speed-insensitive  sources  such  as  oil  seal  and  bearing  preload 
drag  as  well  as  from  speed-sensitive  lubrication  viscous  drag.  Due  to  the 
popular  use  of  hypoid  ring  and  pinion  gear  sets  to  provide  a  lowered  drive- 
shaft  offset,  sliding  gear  viscous  losses  make  the  typical  rear  axle  some¬ 
what  less  efficient  than  a  bevel-gear  set  or  chain  drive  (Ref  42).  Since 
certain  losses  are  relatively  constant  (e.g.  seal  drag),  axle  efficiency, 
defined  by 


P_„  -  losses 

nA  =  __in _  .  ioo 

A  p 

IN 


(  22  ) 


where  P^.^  is  the  input  power  to  the  pinion  shaft,  becomes  a  function  of 
input  power  as  well  as  speed.  Figure  2  shows  an  experimentally  obtained 
typical  axle  efficiency  map  (Ref  1).  Note  that  the  abscissa  coordinate 
is  percent  of  rated  input  power,  indicating  that,  generally,  maximum 
efficiency  occurs  with  minimally  sized  components.  A  least-squares 
quadratic  fit  of  the  data  represented  by  Fig  2  yields  the  following: 

nA  =  93.8  -  .01217 (PSPD)2  +  ,8879(PSPD)  +  4.261  (  23  ) 


where  PSPD  and  PP  are  percent  rated  speed  and  percent  rated  input  power 
respectively.  Eq  (  23  )  is  a  very  accurate  fit  to  the  data  on  Fig  2  with 
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errors  of  less  than  .2  percent  at  all  power  levels  above  15  percent. 


Conventional  fixed-ratio  multi-speed  transmissions  have  efficiency 
maps  very  similar  to  that  of  a  rear  axle  (Fig  3).  Since  hypoid  gears  are 
not  used,  however,  the  sliding  friction  component  is  reduced,  and  consequently 
straight  spur  or  helical-cut  transmission  efficiencies  are  considerably  less 
speed  dependent.  The  efficiency  is  dependent,  however,  on  the  reduction 
ratio  with  1:1  being  the  most  efficient.  Maximum  efficiencies  are  reported 
to  vary  from  99  percent  at  1:1  ratio  to  94.5  percent  at  2.55:1  (Ref  28). 

These  maximum  efficiencies  can  be  closely  approximated  with  the  following 
exponential  equations  as  functions  of  gear  ratio  GR: 

=  99e“,03(GR_1)  (  24  ) 


for  gear  ratios  greater  than  1:1,  or 


nMAX 


99e-.03(l-GR) 


(  25  ) 


for  overdrive  ratios  (i.e.  gear  ratios  less  than  1:1).  As  with  the  rear 
axle  efficiency  derivation,  a  least-squares  quadratic  fit  of  a  published 
transmission  efficiency  map  (Ref  1)  yields  an  accurate  equation  in  terms  of 
percent  rated  shaft  speed  PSPD^,  percent  rated  input  power  PP,  and  maximum 
efficiency  This  equation,  to  be  used  with  Eqs  (  24  )  or  (  25  )  is 

. 00405(PSPD_)2  -  . 1849(PSPD_)  +  1.565 

nT  =  W"  “  - I - ^ - I - )  (  26  ) 


A  variety  of  variable-ratio  transmissions  have  been  suggested  for 
use  with  electric  vehicles  (Refs  2,  3,  19,  32,  36,  37).  These  include 
infinitely-variable  roller  traction  drives,  hydrostatic  drives,  and 
continuously-variable  belt  drives.  Due  to  the  nonlinearities  and 
peculiarities  associated  with  each  of  these  non-standard  transmissions,  no 
attempt  was  made  to  "close-fit"  an  efficiency  equation  for  each  case. 
Instead,  this  model  uses  their  efficiency  data  stored  in  the  form  of  a 
tabular  two-dimensional  performance  map  with  dependent  variables  of  speed 
ratio  and  percent  rated  input  power.  A  look-up  procedure  is  then  used 
with  two-dimensional  linear  interpolation  when  appropriate. 


Motors 


For  reasons  mentioned  in  Chapter  I,  this  vehicle  model  considers 
only  DC  motors.  However,  both  series  and  shunt  motors  are  candidates  and, 
as  such,  are  individually  represented.  Compound  motors  (with  both  shunt 
and  series  field  windings)  are  not  modeled,  but  depending  upon  the  relative 
series  compounding  ratio,  may  be  compared  with  the  behavior  of  a  similarly- 
sized  series  or  shunt  motor.  The  essential  requirement  of  any  motor  model 
is  to  be  able  to  relate  output  power  and  current  for  any  operating  speed 
or  applied  voltage. 

Series  Motor.  In  series  motors,  the  armature  current  I  passes 
entirely  through  the  field  circuit.  Figure  4  shows  a  simplifiei  circuit 
analysis  model. 

The  complete  voltage  equation  for  Figure  4  is 


(  27  ) 


V=1(R.+R)+E  +  (L ,  +  L  )  ~r~ 

M  af  a  a  f  adt 

where  and  are  the  field  and  armature  inductances,  is  the  applied 
motor  voltage,  and  E  is  the  induced  center-electromagnetic  force  (EMF) 
in  the  armature.  The  dynamic  inductance  terms  at  the  end  of  Eq  (  27  ) 
are  significant  only  for  very  short  intervals  of  rapid  current  change  and 
for  this  analysis  can  be  ignored  (Ref  27).  Electromagnet^  power  of  the 
motor  is  defined  as  the  input  electrical  power  minus  all  I  R  copper  losses. 
This  can  be  revealed  by  multiplying  both  sides  of  the  steady-state  version 
of  Eq  (  27  )  by  I  and  rearranging  the  terms 
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=  vMi  -  r(R,  +  R  ) 

M  a  a  f  a 


(  28  ) 


From  the  driving  cycle,  which  determines  transmission  and  controller  position, 
the  motor  speed  n  and  can  be  determined.  Thus  the  variables  on  the 
right  side  of  Eq  (  28  )  must  be  further  reduced  into  functions  of  only 
n  and  to  obtain  an  equation  of  power  useful  in  the  simulation.  Solving 
for  I  ,  again  using  the  steady  state  Eq  (  27  ),  results  in 
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(  29  ) 


R,  +  R 
f  a 


The  counter  EMF  is  proportional  to  the  field  strength  $  and  motor  speed. 
For  series  motors,  assuming  a  linear  magnetization  curve  (i.e.  unsaturated 
fields),  4)  is  proportional  to  I  .  Thus,  from  Eq  (  29  ) 

a  V 

T  M  -  nk<ft  (  30  ) 

a  “  R,+R 
f  a 


VM  -  nk*I 
M _ a 

R,  +  R 
f  a 


(  31  ) 


where  k*  is  a  new  combined  proportionality  constant.  Separating  the 
armature  current  terms  from  Eq  (  31  )  and  simplifying  results  in 


X  =  - - - 

a  Rj  +  R  +  nk* 


(  32  ) 


which  has  the  necessary  and  n  dependence  for  reinsertion  into  Eq  (  28  ). 
The  electromagnetic  power  is  finally  then 


.  Vk* 

(Rf+R  +nk* )' 
t  a 


(  33  ) 


Equations  (  32  )  and  (  33  )  may  be  combined  to  develop  an  overall  motor 

performance  map  as  is  shown  in  Fig  5  for  a  typical  series  motor.  The 
essential  series  motor  characteristics  of  high  torque  (i.e.  power)  at 
low  speed  and  theoretical  unlimited  no  load  speed  can  be  seen  from  this 
figure.  The  actual  operating  envelope  of  any  motor  is  limited  by  the 
applied  voltage  (somewhat  less  than  nominal  battery  voltage,  see  Battery 
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Model,  this  chapter)  on  one  side  and  by  the  maximum  allowable  armature 
current  on  the  other.  Depending  upon  the  robustness  of  the  individual 
series  design,  manufacturers'  recommended  maximum  currents  are  from 
two  to  five  times  the  normal  current  at  rated  power  and  speed  (Refs  4, 
13). 

Shunt  Motors.  The  armature  current  in  a  shunt  field  motor  does 
not  pass  through  the  field  circuit;  in  fact,  for  maximum  flexibility,  the 
field  is  often  separately  excited  as  shown  in  Fig  6.  Neglecting  the 
dynamic  inductance  voltages  as  before,  the  steady  state  circuit  voltage 
equations  are: 


E  +  I  R 
a  a  a 


(  34  ) 


Vf  -  XfRf 


(  33  ) 


where  and  I,  refer  to  the  separate  field  circuit.  Similarly  to  the 
series  motor,  the  electromagnetic  power  equation  may  be  developed  from 
Eqs  (  34  )  and  (  35  ) .  In  this  case,  however,  the  field  flux  <J>  is  pro¬ 
portional  to  the  field  current  which  from  Eq  (  35  )  equals  *  R^.  The 

final  steps  in  the  derivation  result  in 


\  '  ”k*<R7> 
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(  37  ) 


The  shunt  motor  performance  map  on  Fig  7,  developed  from  Eqs  (  36  )  and 
(  37  )  where  =  V  ,  illustrates  the  need  for  separate  field  control  with 
shunt  motors  to  improve  low  speed  performance  and  speed  flexibility.  Again 
the  operating  envelope  is  limited  by  available  voltage  and  maximum  allowable 
armature  current.  Typically,  shunt  armature  currents  should  not  exceed  two 
or  three  times  the  rated  value  to  prevent  commutation  problems  and  over¬ 
heading  (Refs  4,  13). 


Common  Parameters.  The  motor  constants  Ra,  R^,  and  k*  in  Eqs  (  33  ) 
and  (  37)  must  be  resolved  into  functions  of  the  motor's  power  rating  P  to 
provide  the  simulation  with  the  ability  to  rescale  the  motor  to  meet  any 
particular  design  criteria.  The  following  relations  were  developed  for 
modeling  motors  of  different  P  with  the  same  rated  voltage  and  speed: 


U Icct  romag 


where  the  prime  superscript  denotes  changed  values  over  the  basic  P  for 

which  k* ,  r  ,  R  ,  and  I  for  any  operating  condition  are  known  (Refr4l). 
a  i  a 

To  this  point,  all  motor  performance  has  been  based  on  available 
electromagnetic  power,  E  I  .  Additional  mechanical  losses  and  stray-load 
losses  reduce  the  actual  power  output  to  some  fraction  of  the  E  I  (Ref  4). 
These  losses  include  bearing  and  brush  drag,  rotor  windage,  f  iefd^ysteresis 
and  eddy  current-induced  heating.  The  total  effect  is  an  involved  function 
of  motor  speed,  voltage,  and  current  and  is  difficult  to  establish 
analytically.  Available  experimental  data  are  in  the  form  of  overall 
motor  efficiency  curves  which  relate  net  output  power  to  electrical  input 
power  The  motor  model  for  this  simulation  uses  motor  efficiency  data 

nM  stored  in  tabular  form  as  described  in  Appendix  D  to  obtain  net  output 
power  from  the  calculated  EaIfl  by  the  following  equations: 


for  series  motors  or 


for  shunt  motors. 
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(  40  ) 


(  41  ) 


The  final  motor  parameter  required  for  complete  modeling  is  mass 
MM  as  related  to  P  and  rated  speed  n  .  For  constant  voltage  and  speed 
ratings,  the  amounf  of  active  magnetic  material  in  the  field  and  armature 
determines  the  torque  (hence,  power)  of  a  motor.  Reference  27  suggests 
that  a  logarithmic  relationship  generally  exists  between  motor  mass  and 
power  rating.  Indeed,  a  graphical  analysis  of  manufacturer's  data  (Ref  13) 
reveals  the  following  relationship  between  3200  rpm  rated  motors: 


=  {5.75l0s(Pr)}  21.5 


(  42  ) 


A  similar  logarithmic  relation  exists  between  equal  powered  motors  rated 
at  different  speeds.  Again,  analysis  of  data  from  Ref  13  provides  the 
following  version  to  Eq  (  42  )  for  motors  rated  at  other  than  3200  rpm: 


=  {5.75  1°8(Pr)}21.5{.162log(37D0)}  (  43  ) 


Controllers 


Motor  controllers  must  regulate  the  motor  voltage  to  adjust  power 
output  and  limit  armature  current  at  low  speeds.  To  be  satisfactory,  they 
must  accomplish  this  control  with  minimum  power  loss. 


Step-Voltage  Switching.  Motor  controllers  which  employ  mechanical 
relays  to  adjust  the  parallel-series  battery  network  dissipate  power  only 
through  the  relay  activator  circuits  and  slight  contact  resistance.  Properly 
designed,  their  efficiencies  can  exceed  99  percent  (Ref  32).  Ignoring  this 
negligible  power  loss,  modeling  consists  of  simply  knowing  the  battery 
voltage  fraction  available  at.  each  step  and  the  corresponding  parallel- 
series  network.  Analysis  of  the  net  source  impedance  requires  this  know- 


ledge  of  the  specific  battery  arrangement  for  each  step.  Since  the  number 
of  individual  relays  required  increases  dramatically  with  the  number  of 
steps  provided  (e.g.  nine  relays  for  three  steps,  21  relays  for  four  steps, 
(Ref  30),  the  practical  upper  limit  for  these  controllers  is  three. 

To  accommodate  additional  voltage  steps  while  reducing  the  relay 
requirement,  high  current  diodes  may  be  inserted  to  replace  selected  relays. 
The  simplicity  and  economy  they  offer  over  relays  is  partially  balanced  by 
the  losses  they  create.  Depending  upon  the  semi-conductor  material  used  in 
their  construction,  these  diodes  will  typically  produce  about  a  one-volt 
drop  when  normally  conducting  with  forward  bias.  Discriminate  placing  of 
these  diodes  minimizes  their  dissipative  losses  by  requiring  only  a 
fraction  of  the  total  current  to  pass  through  any  single  diode.  The  net 
effect  on  the  switching  controller  model  by  the  addition  of  diodes  is  to 
lower  each  step's  voltage  fraction  by  the  appropriate  diode-induced  voltage 
drop. 

Solid-State  Inverters  (Choppers).  Electronic  chopper  controls  offer 
continuous  voltage  control  from  zero  to  95  percent  of  the  battery  voltage. 
Their  basic  principle  is  to  switch  the  motor  on  and  off  at  a  high  frequency 
(up  to  1000  Hz)  while  controlling  the  relative  duration  of  the  conduction  and 
nonconduction  periods  to  produce  varying  time-averaged  voltages.  The  power 
losses  associated  with  choppers  result  from  the  resistive  and  semi-conductor 
elements  used  in  the  oscillator  circuits.  The  combined  losses  have  been 
modeled  (Ref  2)  as 

P.  *  (1.5  volts ) +  ( .03)E  I  (  44  ) 

C  M  a  a 

where  the  first  term  on  the  right  side  is  the  semi-conductor  junction  loss 

and  the  second  term,  the  resistive  losses  as  a  function  of  motor  power. 

These  losses  can  be  more  easily  accounted  for  when  referred  back  to  a 

voltage  drop.  Dividing  Eq  (  44)  by  the  supply  current,  assuming  that 

Iw  *  I  and  that  E  =  V„,  results  in 
M  a  a  M 


Vc  -  1.5  +  .03Vm 


(  45  ) 


This  controller  loss  then  reduces  the  net  instantaneous  battery  voltage 

(  46  ) 


V  available  to  the  motor  since 

u 


V  =  V  -  V 
M  B  C 


or  combining  Eqs  (  45  )  and  (  46  ) 


V  =  VB  "  X-5 
M 


(  47  ) 
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Batteries 


While  numerous  battery  types  promising  badly-needed  energy  and  power 
density  improvements  lead-acid  batteries  with  their  commercial  availability 
and  proven  reliability  form  the  basis  for  this  battery  model.  The  model 
could  be  modified  for  alternate  battery  types  by  substituting  the  appropriate 
performance  parameters.  Currently,  however,  insufficient  published 
experimental  data  are  available  to  detail  the  dynamic  performance  of  these 


batteries,  although  the  Shepherd  Equation  (Refs  38,  39),  with  its 
empirically-derived  constants,  provides  some  insight.  Any  useful  battery 
model  must  be  able  to  predict  power  available  and  energy  remaining  at 
any  time  during  a  driving  cycle  as  functions  as  battery  mass,  ambient 
temperature,  discharge  level,  load  current,  nominal  voltage,  and  loading 
history. 


Energy  Capacity.  The  constant  load  useful  capacity  of  a  battery 
is  a  function  primarily  of  discharge  rate.  Figure  8  shows  the  constant 
load  capacity  characteristics  (Ref  32)  used  in  this  model.  A  least-squares 
quadratic  approximation  of  that  data  results  in 

c  =  195  -  . 633 C Ig )  +  9.13  •  10-4(Ig)2  (  48  ) 

where  battery  capacity  c  is  expressed  in  amp-hrs  and  I  is  the  steady 
battery  discharge  current.  The  capacity  of  lead-acid  batteries  is  also 
a  strong  function  of  temperature  (Ref  43).  Eq  (  48  )  is  correct  only  for 
the  standard  test  temperature  of  25  C.  For  other  temperatures,  the  capacity 
is  modified  by  factor  c* ,  where 


c*  =  .014 (T)  +  .65  (  49  ) 

Incorporating  this  into  Eq  (  50  )  results  in 

c  =  {195  -  .633(Ib)  +  9.13  *  10"4(Ig)2} 

•  { .014(T)  +  .65}  (  50  ) 

Eq  (  50  )  is  still  only  correct  for  constant  IR  values.  An  incremental 
scheme  based  on  sequentially  solving  this  equation  for  each  incremental 
value  of  Ig  would  yield  an  incorrect  total  discharge  value. 

During  discharge  of  any  battery  with  liquid  electrolyte,  ion 
migration  toward  their  respective  electrodes  within  each  cell  causes  the 
output  to  be  time-dependent  (Ref  38).  This  developed  non-homogeneity, 
in  effect,  temporarily  reduces  the  capacity  of  the  battery.  Under  lighter 
loads  or  no  load,  the  electrolyte  eventually  redistributes  itself  homo¬ 
geneously  again.  The  consequence  of  this  time-dependent  characteristic  is 
that  the  instantaneous  capacity  of  a  battery  is  a  strong  function  of 
previous  discharge  rate  as  well  as  present  load.  Likewise,  idle  time 
actually  provides  a  recuperative  effect  to  marginally  regain  capacity.  The 
following  discharge  equations  consider  load  history  in  terms  of  average 
current  Ig  and  have  shown  good  correlation  with  experimental  results  (Ref  6): 

DF  =  {~B  *  C  (  B)  +  XB  '  C(1  -  IB)}3600  (  51  ) 


for  instantaneous  current  ID  <  average  current  I  or 

O  D 
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(  52  ) 


DF  =  {IB  *  t(IB)  +  IB  ‘  C  (1  -  IB)}3600 


for  Ig  >  I  .  The  discharge  fraction  DF  represents  the  accumulated  portion 
of  the  battery  used  during  the  cycle  time  t  while  c  and  T  are  the  solutions 
to  Eq  (  50  )  based  on  IR  and  IR  respectively.  Since  Eq  (  50  )  was  developed 
from  data  for  a  single  standard  size  (~  29  kg),  six-volt,  golf-cart  style 
battery,  any  simulated  change  in  battery  mass  must  appropriately  affect  the 
discharge  life.  This  model  assumes  that  increased  battery  mass  results  in 
the  addition  of  electrode  plate  area,  hence  the  current  IR  in  Eqs  (  51  )  and 
(  52  )  may  be  replaced  by  a  current  I  '  such  that 

s’  •  < 53  > 

where  BM'  is  the  new  battery  mass  and  BM  is  the  nominally-sized  mass  based 
on  the  batteries  characterized  by  data  given  in  Fig  9. 

Available  Power.  The  output  of  a  battery  is  equal  to  its  output 
voltage  times  the  load  current.  The  load  current  is  a  function  of  the 
vehicle  demands  while  the  output  voltage  depends  upon  that  current  as  well 
as  the  instantaneous  discharge  state  of  the  battery  and  electrolyte 
temperature.  For  constant  temperature,  the  battery's  volt-amp  characteristic 
may  be  modeled  as  a  constant  resistance  (Ref  1,  32,  38,  43).  For  the  lead- 
acid  type  considered,  this  resistance  is  assumed  to  be  .003  ohms  per 
standard  six-volt  battery  at  25  C.  As  with  capacity,  this  resistance  value 
is  also  a  function  of  temperature  (Ref  43).  A  quadratic  approximation  of 
this  temperature  dependence  produces  the  following  equation  for  relative 
resistance 

Rg*  =  ( .00087 )T2  -  (.047)T  +  1.63  (  54  ) 

Hence  the  overall  resistance  per  battery  is 

Rg  =  ( .003 )Rg*  (  55  ) 

Battery  connector  cable  impedance  is  not  explicitly  modeled  since  properly 
sized  wire  with  short  lengths  will  keep  this  value  an  order  of  magnitude 
below  the  battery  resistance.  However,  when  significant,  this  resistance 
may  be  added  to  Eq  (  49). 

A  final  consideration  in  the  voltage  profile  of  a  battery  is  the 
present  stage  of  discharge.  In  general,  the  voltage  produced  at  any 
temperature  and  load  will  decrease  as  the  battery  becomes  increasingly 
discharged.  For  this  model,  in  fact,  100  percent  state-of-discharge  is 
defined  as  when  the  cell  potential  has  dropped  to  1.75  volts  at  a  steady 
load  of  80  amps.  This  definition,  combined  with  typical  lead-acid 
discharge  characteristics  (Refs  1,  14),  results  in  the  following  equation 
for  no-load  cell  voltage: 

V  =  2.03  +  ( .08)DF  -  (.28)DF2  (  56  ) 

c 
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Including  resistance  effects  (Eq  49),  the  general  voltage  expression  for 
a  nominal  six-volt  battery  thus  becomes: 

V„  =  6.09  -  ( . 003)R*I„  +  ( . 24 )DF  -  (,84)DF2  (  57  ) 

B  D  D 

To  model  a  change  in  battery  mass,  the  nominal  resistance  Rg  is  replaced  by 
a  new  value  R^'  related  to  the  BM  change  by 

^  =  ^  58  ^ 

where  BM'  is  again  the  new  battery  mass. 

An  added  consideration  for  a  useful  battery  model  is  the  net  battery 
efficiency.  This  is  defined  as  the  ratio  of  energy  discharge  to  the  charging 
energy  (Ref  22).  Therefore,  this  efficiency  is  a  function  not  only  of  dis¬ 
charge  rate  and  extent,  but  also  the  charging  rate  and  extent  of  charge. 

In  general,  high  discharge  and  charge  rates  lower  the  efficiency  due  to 
internal  resistance  losses  and  hydrogen  gas  formation.  However,  the  most 
efficient  rate  of  charge  varies  directly  with  the  battery's  depth  of  charge. 
Also,  due  to  internal  polarization,  some  experimental  evidence  indicates 
that  high  charging  rates  immediately  following  high  discharge  rates  can 
be  very  efficient  (Ref  1).  In  view  of  the  complexity  of  dynamically  modeling 
both  charge  and  discharge  cycles,  this  battery  model  uses  an  experimentally 
derived  (Ref  32)  average  efficiency  of  70  percent. 

Hybrid  Power.  As  discussed  in  Chapter  I,  the  major  utility  of  gas- 
powered  hybrid  augmentation  is  to  relieve  battery  drain  during  high  speed 
cruise.  For  that  reason,  the  hybrid  is  modeled  as  being  inactive  until  a 
designated  engagement  speed.  At  that  speed  the  hybrid  becomes  positively 

connected  to  the  motor  shaft  via  some  transmission  arrangement  and  is 

capable  of  producing  any  power  up  to  the  maximum  for  its  particular  operating 
speed.  The  hybrid  engine  currently  included  in  the  vehicle  model  is  a  7.5 
kw,  two-cycle  industrial  engine  (Ref  45).  From  published  performance  curves, 
the  maximum  power  available  from  this  hybrid  engine  is: 

PR  =  {(8.7  •  10"4)nH  +  .99}nR  (  59  ) 

where  the  engine  speed  equals  the  motor  speed  multiplied  by  the  hybrid 

transmission  ratio  and  is  the  average  power  efficiency  of  that  trans¬ 
mission.  The  fuel  consumption  of  the  engine  is  also  taken  from  manu¬ 
facturer's  data  (Ref  23)  and  is  modeled  as  a  constant  brake-specific  fuel 
consumption  of  .56  liters  per  kwh. 

This  final  entry  now  completes  the  vehicle  model  characteristics 
required  to  obtain  the  net  input  energy  necessary  to  meet  the  driving  cycle 
road  load  power  demands.  It  remains  for  the  simulation  then  to  combine  all 
the  dynamic  characteristics  of  a  particular  vehicle  system  with  a  driving 
and  speed  schedule  to  derive  the  instantaneous  and  overall  vehicle 
performance . 
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III.  Simulation  Methodolog 


Figure  9  is  a  simplified  flow  chart  of  the  simulation  program  EVSIM 
developed  during  this  effort.  Reference  41  includes  a  user's  guide  and 
includes  a  sample  computer  listing,  data  deck,  and  sample  output.  Additional 
information  on  two  of  the  more  involved  subroutines,  VSTF.P  and  TRANSEFF,  are 
also  located  in  Reference  41. 

The  simulation  begins  by  selecting  certain  "fixed"  vehicle  and 
environment  parameters.  These  represent  a  particular  vehicle  type  and  sur¬ 
rounding  influences  and  include  such  things  as  payload  mass,  transmission 
type,  and  motor  type  as  well  as  ambient  temperature  and  fuel  costs.  Certain 
"basic"  vehicle  parameters  are  also  entered  at  the  start.  These  are  parame¬ 
ters  such  as  motor  power  and  battery  mass  which  may  be  incrementally  changed 
by  the  simulation  to  meet  specific  performance  criteria.  Table  II  lists  all 
the  specific  input  parameters  grouped  by  their  category. 

An  input  driving  cycle  is  needed  to  provide  a  second-by-second  speed 
schedule.  For  each  time  increment,  road  load  forces  are  calculated  and  con¬ 
verted  to  power  required  at  the  wheels.  Through  continuous  efficiency 
equations  and/or  component  efficiency  maps,  this  raw  road  power  is  converted 
to  the  required  motor  output  power.  Calculation  of  the  available  motor  power 
depends  upon  the  controller  type.  For  step-type  controllers,  only  discreet 
power  steps  are  available  for  any  particular  motor  speed  (Reference  41); 
hence,  an  exact  match  with  the  required  power  at  any  increment  is  unlikely. 
The  power  step  closest  to  the  power  required  is  selected  and  then  an  iter¬ 
ative  procedure  accomplished  by  adjusting  the  interval's  terminal  speed  until 
the  power  required  does  match  the  available  power.  For  continuous  con¬ 
trollers,  EVSIM  assumes  a  perfect  match  between  required  and  available  power 
until  the  maximum  available  power  is  insufficient  to  meet  the  demands.  Then, 
as  with  the  step-controller  an  iterative  procedure  determines  the  actual 
speed  attained  at  the  end  of  each  interval.  For  systems  using  a  gas-engined 
parallel  hybrid,  the  simulation  assumes  that  the  continuously  controllable 
hybrid  power  is  used  only  to  match  the  total  available  power  with  that 
required . 


Once  the  actual  vehicle  speed  and  motor  power  for  each  time  increment 
have  been  resolved,  the  fraction  of  battery  discharge  is  calculated.  This 
entire  process  is  repeated  until  the  cycle  is  complete  or  until  the  vehicle 
fails  to  meet  some  minimum  critical  performance  criterion.  If  the  basic 
vehicle  is  insufficiently  powered,  an  excessive  speed  deviation  from  the 
driving  cycle  schedule  will  cause  EVSIM  to  increase  the  rated  power  and 
associated  parameters  until  the  schedule  can  be  achieved.  On  the  other  hand, 
if  lack  of  energy  storage  is  apparent  from  an  excessive  discharge  fraction, 
EVSIM  can  increase  the  battery  mass.  Naturally,  both  of  these  parameter 
adjustments  are  interrelated  and  typically  one  will  require  the  other. 

When  the  input  driving  cycle  has  been  successfully  completed  the 
final  outputs  of  EVSIM  include  the  total  motor  energy  used,  the  total  deceler¬ 
ation  energy  available  at  the  motor  (for  regeneration),  the  percent  of  battery 
discharge  at  cycle  completion,  the  average  electrical  cost  per  kilometer,  and 
the  final  vehicle  parameters  if  adjustments  were  made.  The  power  cost  is 
based  on  an  average  charger  efficiency  of  85  percent  (Reference  32)  combined 
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with  the  battery  model's  70  percent  efficiency  to  produce  a  net  60  percent 
energy  transfer  from  AC  power  at  the  charger  to  DC  power  out  of  the  batteries 
If  a  parallel  hybrid  system  is  simulated,  EVSIM  also  recalls  the  total  engine 
energy  used  along  with  the  average  petroleum  fuel  cost  per  kilometer. 


TABLE  II 


EVSIM  Input  Parameters 


FIXED 

BASIC 

Drag  coefficient 

Payload  mass 

Fixed  mass 

Frontal  area 

Rolling  radius 

Rolling  coefficient 

Axle  ratio 

Transmission  type 

(Multi-ratio) 

gear  ratios 
shift  speeds 

(CVT ) 

efficiency  data 

_  , 

Motor  type 

Rated  power 

base  speed 

Armature  resistance 

efficiency  data 

Field  resistance 

maximum  power  ratio 

k*  value 

Maximum  current 

Controller  type 

Step  source  impedance 

(Step  voltage) 

voltage  fractions 

Nominal  battery 
voltage 

Hybrid 

Battery  mass 

coupling  efficiency 
speed  ratio 
engagement  speed 
operating  gear 

Ambient  temperature 

Electric  energy  cost 
Petroleum  fuel  cost 

Speed  schedule 

Program  option 

A  critical  feature  in  any  vehicle  simulation  is  the  realism  of  the 
driving  sequence.  Several  test  cycles  have  been  devised  for  both  con¬ 
ventional  and  electric  vehicles.  While  the  more  elaborate  cycles  more  accu¬ 
rately  represent  actual  driving  experiences,  the  tradeoff  is  the  increased 
difficulty  in  simulating  these  sequences  and,  consequently,  obtaining  vehicle 
road  loads . 


fcle  Selection 


In  February  1976  the  SAE  J227a  Electric  Vehicle  Test  Procedure  was 
developed.  This  cycle  was  designed  primarily  for  track  testing  and  was  pro¬ 
duced  in  three  versions,  the  most  demanding  of  which  is  shown  in  Figure  ID. 
This  schedule,  with  its  top  speed  of  72  km/hr,  is  most  frequently  used  to 
simulate  "variable-route  suburban"  driving.  As  only  the  end  points  of  each 
portion  are  specified,  distinct  acceleration  or  d“~  .deration  profiles  are  not 
required.  The  relative  ease  of  simulating  this  driving  cycle  has  led  to  its 
frequent  use  (References  1,  9,  19,  32,  37).  However,  when  used  in  track 
tests  (Reference  32),  it  generally  overpredicts  range  and  underestimates 
energy  usage  rate  when  compared  to  user  experience. 


A  much  more  intricate  driving  cycle  is  the  Federal  Urban  Driving 
Sequence  (FUDS,  Reference  10),  which  was  developed  primarily  for  dynamometer 
vehicle  emission  tests.  It  is,  however,  also  used  by  the  Environment  Pro¬ 
tection  Agency  for  new  car  urban  fuel  economy  tests.  This  23-minute  cycle  is 
the  recorded  operation  of  an  actual  vehicle  being  driven  in  the  Los  Angeles 
area.  As  such,  it  is  a  second-by-second  speed  schedule  that  includes 
frequent  stops,  random  traffic-induced  speed  adjustments,  and  one  short 
stretch  of  90  km/hr  freeway  travel  (Figure  11).  In  spite  of  its  complexity 
(1371  data  points),  its  realism  and  ability  to  provide  a  base  for  operating 
cost  comparison  with  new  conventional  vehicles  made  the  FUDS  the  logical 
choice  for  this  simulation. 


Simulation  Output 


The  sample  output  on  Figure  12  demonstrates  EVSIM's  evaluation  mode. 
The  vehicle  modeled  has  the  configuration  of  the  test-bed  vehicle  design.  It 
represents  a  series  motor/gasoline  engine  parallel  hybrid.  The  transmission 
is  a  four-speed  conventional  gear  ratio  model  (although  only  the  first  three 
gears  are  used  during  hybird  augmentation).  The  controller  is  a  five-step 
battery-switching  system  that  uses  a  combination  of  relays  and  high  power 
diodes  to  control  the  motor  supply  voltage.  The  selected  portions  of  the 
complete  cycle  included  in  Figure  12  indicate  the  extent  of  computed  data 
available  for  each  increment  of  the  chosen  driving  cycle.  The  final  page  of 
sample  output  shows  the  overall  energy  recap  that  is  computed  at  the  end  of 
each  completed  simulation  cycle. 


EVSIM's  cruise  sequence  at  a  constant  88  km/hr  is  demonstrated  in 
Figure  13.  The  same  performance  and  energy  use  summary  is  presented 
following  the  complete  discharge  of  the  batteries. 
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Federal  Urban  Driving  Sequence 


•  ♦♦ELECTRIC  VEHICLE  SI"ULATION*** 


VEHICLE  CONSTANTS  AM: 

CRAG  COEFFICIENT-  .A2 
FRCMAL  AREA-  1.B0  SCUARE  PETERS 
PAYLCA0  PASS-270.  K 1LOCRAPS 
FIXEC  PASS-325. 

ENVIRCN*ENTAL  CONSTANTS  ARE: 

APBIENT  TEPPERATURE-  27.  CELSIUS 
ELECTRIC  ENERGY  COST-  S.C50  PER  KWHR 
PETRCLEUP  FUEL  COST*  1  .35  PER  LITER 

POWER  REDUCTION  VARIABLES  ARE: 

TIRE  ROLLING  RADIUS-  -  3  *•  3  PETERS 
ROLLING  COEFFICIENT-  .0105 
AXLE  RATIO-  A. 13 
NUPPER  FORWARD  gears-  4 
PAXIPIP  VELOCITY-112.  KP/HR 

FIXED  SPEED  TRANSMISSION  VARIABLES  ARE: 

PCTOR  UPSHIFT  SPEEC-1500.  RPP 
w  C  T  C  F  0CWNSH1FT  SPEEC-  500.  RPP 
GEAR  RATIOS*  3.800 
__  2.060 

1.260 

.880 


■C T C R  VARIABLES  ARf: 

POT  OR  TYPE-  SERIES 

RATEC  POWER-  6.0  *W 

PASE  SPEED-3200.  RPP 

PAX  ARPATURE  CURRENT-  400.  APPS 


IPPECANCE-  .  CO  30  OHPS 
.COA 
.COT 
.018 
•  .  C  36 


BATTERY  VARIABLES  ARE: 

NOMINAL  VOLTAGE-  72. 
eATTERY  PASS-  350. 


5  STEP  CONTROLLER  VARIABLES  ARE: 

VOLTACE  FRACTION-  .167  SOURCE 
.333 
•  SCO 
.667 
l.COO 


7KW  SERIES  HYBRID  VARIABLES  ARE: 

HYPR  IC/"CTOR  SPEED  RATIO-  1.950 
COUPLING  EFFICIENCY-  96.5* 
ENCACEPENT  SPEEO-  AS.OKP/HR 
TRANSPISSICN  PDSITICN-  3  GEAR 


SY5TE"  EVALUATICN  OPTION  SELECTEO 


Fig  12  EVSIM  Sample  Output 
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Hybrid  Vehicle  Cruise  Simulation 


Hardware  Implementation  of  the  EV 


Construction  of  a  test  bed  EV  was  needed  as  a  source  of  data  to  vali¬ 
date  the  simulation  EVSIM  and  to  investigate  new  efficiency  and  operability 
concepts  for  EV's.  Since  portions  of  the  preliminary  design  of  the  test  bed 
vehicle  proceeded  in  parallel  with  the  simulation  development,  some  com¬ 
ponents  selection  was  based  on  a  preliminary  analysis  of  the  desired  per¬ 
formance  and  candidate  drive  systems. 

Candidate  Selection 

Table  III  lists  the  major  components  of  any  electric  vehicle  system 
and  includes  the  options  considered  for  the  test  vehicle.  In  addition  to 
the  performance  requirements  listed  in  the  previous  section,  this  vehicle  was 
also 


TABLE  III 

Electric  Vehicle  Components 


Batteries  Controllers  Motor  Transmissions 


Lead-Acid 
-Std  Automotive 
-Industrial 
-Golf-cart 
Nickel  Zinc 
Nickel  Iron 


Battery  Switching  DC  Shunt 

Chopper  DC  Series 


Fixed  Ratio 

Gear  Ratio 

Continuously 

Variable 


constrained  to  use  readily  available  components  and  remain  flexible  enough  to 
incorporate  various  modifications  during  follow-on  investigation.  Many  con¬ 
ceivable  component  options  were  not  included  in  Table  III  due  to  their 
current  lack  of  commercial  availability  or  apparent  lack  of  suitability  for 
the  test  vehicle.  Alternative  current  (AC)  motors  were  not  considered,  for 
example,  due  to  the  control  complications  and  associated  losses  with  the 
required  DC  to  AC  inverters  (Reference  16).  Similarly,  several  highly 
regarded  battery  types  (such  as  zinc-chlorine  hydrate  and  sodium  sulfur), 
although  possessing  significantly  better  potential  performance  than  lead-acid 
batteries  (Reference  14 ,  22,  32,  35),  were  not  options  for  the  test  vehicle 
due  to  their  present  lack  of  development  and  availability. 

While  in  general,  the  interactions  between  the  various  basic  power 
train  components  encourages  a  systems'  viewpoint  toward  their  selection  or 
elimination,  certain  individual  component  decisions  were  made  based  simply 
on  the  constraints  of  the  test-bed  vehicle.  The  reasonable  flexibility  and 
the  immediate  availability  of  a  multi-speed  gear  ratio  transmission  made  it 
the  obvious  design  choice.  The  battery  selection  also  could  be  reasonably 
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separated  from  other  component  influences.  The  21  kwh  energy  storage  and  22 
kw  peak  power  demands  were  both  severe  restrictions  that  eliminated  certain 
battery  candidates. 


With  lead-acid  batteries,  high  energy  and  high  power  densities  cannot 
be  simultaneously  optimized  within  the  same  design  (References  25,  35). 
Standard  automotive  starting,  lighting  and  ignition  (SLI)  batteries,  for 
example,  have  very  high  power  densities  (up  to  200  w/kg  for  starting  inter¬ 
vals,  Reference  32),  but  were  eliminated  from  final  consideration  due  to 
their  inability  to  recover  from  repeated  deep  discharges.  At  the  other  end 
of  the  spectrum  are  the  industrial  designs.  These  are  designed  for  appli¬ 
cations  such  as  electric  lift  trucks,  where  weight  and  physical  size  are  not 
critical  considerations  and  maximum  peak  power  demands  are  not  appreciably 
above  average  levels.  They  can  sustain  up  to  2000  deep  cycle  discharges  that 
have  limiting  power  densities  below  20  w/kg  (Reference  32).  Golf-cart 
batteries  were  developed  to  provide  reasonable  power  densities  (up  to  100 
w/kg)  while  maintaining  deep-discharge  life  expectancies  of  200  to  400 
cycles  (References  22,  32);  as  such,  they  represent  the  best  lead-acid 
alternative . 

Both  high  power  density  and  high  energy  density  are  potentially 
attainable  with  the  two  nickel-based  battery  candidates.  Unfortunately, 
recent  laboratory  and  vehicle  test  of  nickel-zinc  designs  have  failed  to  sub¬ 
stantiate  their  projected  performance  and  have  shown  especially  poor  deep- 
discharge  cycle  life  (Reference  15).  The  nickel-iron,  or  "Edison"  battery, 
on  the  other  hand,  has  demonstrated  excellent  cycle  life  (1000+),  as  well  as 
high  power  density  (130  w/kg)  and  approximately  25  percent  higher  energy 
density  than  lead-acid  batteries.  Poor  energy  efficiency,  however,  is  a 
strong  deterrent  against  the  design.  Due  largely  to  excessive  hydrogen  gas 
evolution  during  charging,  average  charge/discharge  efficiency  is  typically 
50  percent  which  is  about  one-third  lower  than  that  reported  for  lead-acid 
batteries.  Considering  all  the  reasonable  options  then,  lead-acid  golf-cart 
style  batteries  were  selected  as  being  the  most  practical  and  energy 
efficient . 
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The  choice  of  a  motor  and  controller  was  not  as  clear.  Series  motors 
have  traditionally  been  used  for  traction  motors  due  to  their  high  starting 
torque  characteristic.  Shunt  motors,  on  the  other  hand,  offer  much  better 
speed  regulation  and  the  opportunity,  with  separate  field  control,  to  greatly 
reduce  the  current  handling  requirements  of  the  control  circuitry.  Either  of 
the  proposed  controller  types  can  be  effective  with  either  motor;  however, 
the  over-current  sensitivity  of  the  shunt  motors  at  low  speeds  (Reference  16) 
strongly  encourages  the  sort  of  current-limiting  circuits  as  are  typically 
included  in  solid  state  chopper  controls.  In  the  final  analysis,  the  low 
speed  torque  capability  of  the  series  motor  and  the  simplicity  and  flexi¬ 
bility  of  the  battery-switching  controller  prompted  their  selection  for 
incorporation  into  the  test-bed  vehicle. 


Test  Vehicle  Design 


Once  the  general  configuration  of  lead-acid  batteries,  battery 
switching  controller,  series  motor,  and  multi-ratio  transmission  was 


selected,  the  actual  component  sizing  and  vehicle  designing  could  proceed. 

The  most  obvious  physical  constraints  came  directly  from  the  vehicle  chassis 
that  was  available.  The  total  number  of  batteries  that  could  be  conveniently 
installed  within  the  body  of  the  car  while  retaining  four  passenger  capa¬ 
bility  was  twelve.  Based  on  the  standard  29  kg,  six-volt  golf-cart  battery, 
this  constraint  resulted  in  a  72  volt,  350  kg  battery  pack.  The  net  energy 
capacity  of  the  system  at  a  two-hour  discharge  rate  was  calculated  by  multi¬ 
plying  the  energy  density  of  30  wh/kg  (Reference  32)  by  the  battery  mass  to 
obtain  10.5  kwh.  This  realization  pointed  out  the  facts  that  the  battery 
system  alone  could  not  provide  the  21  kwh  necessary  for  the  desired  160  km 
cruise  range,  and  that  the  maximum  power  that  could  be  drawn  from  the  batter¬ 
ies  over  the  two  hours  was  5.25  kw.  This  power  capacity  led  to  motor  power 
rating  of  6  kw.  Since  series  motors  can  develop  up  to  four  times  their  rated 
power  (References  4,  13)  for  limited  durations,  the  22  kw  acceleration  power 
requirements  was  attainable. 

The  most  significant  conclusion  of  the  battery  energy  limitation  was 
that  some  auxiliary  energy  source  (i.e.  hybrid)  was  required.  Two  gasoline- 
powered  hybrid  arrangements  were  considered.  A  series  hybrid  system  was 
first  investigated,  but  disgarded  due  to  its  low  conversion  efficiency  of 
approximately  forty  percent  of  mechanical  motion  to  stored  battery  energy 
(References  32,  41). 

A  parallel  hybrid  system,  as  indicated  in  Figure  14  reduces  consider¬ 
ably  the  total  losses  to  the  engine/motor  coupling  inefficiencies.  Efficien¬ 
cies  of  approximately  ninety  percent  are  achieveable.  This  was  the  system 
selected  for  the  test  vehicle.  Since  the  motor  could  provide  enough  power  to 
meet  the  urban  acceleration  demands,  the  prime  function  of  the  hybrid  engine 
was  simply  to  augment  the  electric  drive  during  high  speed  cruise.  Besides 
supplying  the  added  power  required  for  curise  speed,  with  appropriate  engine/ 
motor  ratios,  hybrid  augmentation  could  actually  reduce  the  motor  current 
demand,  thus  extending  the  battery  range.  A  coupling  ratio  was  thus 
established  that  would  allow  high  speed  cruise  in  a  relatively  low  trans¬ 
mission  ratio  to  cause  a  high  motor  shaft  speed.  Increasing  the  speed  of  a 
series  motor  automatically  reduces  its  current  demand. 

For  the  test  vehicle,  a  parallel  shaft  V-belt  drive  was  designed. 
V-belt  drive  had  the  advantages  of  high  efficiency  (~  96  percent)  and  speed 
ratio  flexibility  through  alternate  sheave  combinations  (Reference  20).  One 
sheave  was  mounted  onto  the  motor  shaft  such  that  the  hybrid  power  could  be 
directly  coupled  to  the  motor  output,  A  belt-tightener  clutching  arrangement 
permitted  the  hybrid  engine  to  be  disconnected  during  slow-speed  pure 
electric  operating  modes.  The  constantly  turning  motor  sheave  could  then  be 
used  to  start  the  hybrid  engine  when  the  appropriate  vehicle  speed  was 
reached.  For  initial  test  purposes,  the  clutching  mechanism  was  designed  to 
be  manually  operated  by  the  vehicle  driver. 

Once  the  preliminary  design  concepts  were  established  for  the  test 
vehicle  from  use  of  EVSIM  and  analysis,  the  detailed  component  hardware  was 
designed.  Included  in  this  design  was  a  lightweight  motor  mounting  bracket 
that  provides  positive  shaft  alignment  as  well  as  a  mounting  pad  for  the 
hybrid  engine  bracket  and  one  for  any  future  ancillary  test  equipment.  The 
motor  to  transmission  coupling  was  designed  with  an  outside  taper  such  that  a 
variety  of  stock  taper-lock  V-belt  sheaves  can  be  mounted.  This  coupling 


also  incorporates  a  separate  bearing  to  absorb  the  V-belt  created  side  loads. 
A  final  notable  design  feature  for  the  baseline  test  vehicle  was  the  five- 
step  battery-switching  controller  developed.  Using  eight  high-current  diodes 
and  seven  relay  switches,  five  voltage  ranges,  including  12,  24,  36,  48,  and 
72  volts,  are  available. 

Test  Vehicle  Baseline  Hardware  Requirements 

The  performance  requirements  for  the  test-bed  vehicle  were  90  km/hr 
top  speed,  0  to  50  km/hr  in  fifteen  seconds,  two  hours  urban  endurance,  and 
160  km  range  at  80  km/hr  cruise.  The  selected  baseline  system  as  described 
by  EVSIM  in  Figure  12  whose  urban  performance  is  also  indicated  in  Figure  12, 
meets  these  requirements.  Note  from  the  input  parameter  listing  in  Figure  12 
and  the  portions  of  the  driving  cycle  displayed  that  the  selected  hybrid 
system  operates  in  a  pure  electric  mode  until  the  gasoline  engine  engaged  at 
45  km/hr  to  improve  the  high-speed  cruise  performance.  With  the  system's  6 
kw  series  motor,  7.5  kw  gasoline  engine  and  72  volt,  350  kg  lead-acid  battery 
pack,  EVSIM  predicted  a  top  speed  of  90  km/hr  (hybrid  engine  rpm-limited)  and 
zero  to  50  km/hr  acceleration  in  13  seconds.  A  constant  88  km/hr  cruise 
sequence  (Figure  13)  indicated  a  range  of  173  km  for  100  percent  battery  dis¬ 
charge  with  an  ending  top  speed  of  86.9  km/hr.  The  extrapolated  urban 
endurance  based  on  the  battery  discharge  during  one  23-minute  cycle  was  two 
hours  and  37  minutes.  Average  direct  operating  energy  costs  predicted  were 
$. 012/km  for  the  urban  environment  and  $. 0188/km  for  the  high-speed  cruise 
sequence.  This  is  based  on  a  cost  of  $. 35/liter  for  gasoline  and  $.06/kwhr 
for  electricity.  This  differential  can  be  explained  by  the  higher  cost  of 
fuel  being  consumed  by  the  hybrid  engine  during  cruise  as  opposed  to  the 
nearly  pure  electric  power  used  in  urban  operation. 

EVSIM's  system  design  mode  was  also  used  to  develop  an  all-electric 
series  motor  system  that  would  have  the  equivalent  high-speed  range  and  per¬ 
formance  of  the  hybrid  design.  This  configuration  required  1200  kg  of 
battery  and  a  15  kw  motor  to  attain  the  160  km  range.  Additionally,  the 
urban  direct  operating  costs  increased  to  $. 0188/km  while  the  cruise  costs 
marginally  decreased  to  $. 0151/km. 

Once  a  baseline  test  bed  vehicle  was  designed,  it  was  implemented  in 
hardware.  It  was  instrumented  to  monitor  energy  used  (instantaneous  voltage 
and  current,  and  accumulated  watts  used)  speed,  electric  motor  rpm,  and 
internal  combustion  engine  (hybrid)  rpm.  With  this  instrumentation,  the 
driving  cycles  used  in  EVSIM  could  be  easily  validated. 


Validation 

The  EVSIM  simulation  was  validated  from  primarily  two  sources;  the 
test  bed  vehicle  and  available  data  (References  1,  9,  31,  32,  34).  The  availa 
ble  data  accumulated  by  DOE  and  NASA  was  useful  in  verifing  and  validating 
certain  subsystem  models  such  as  the  vehicle  retarding  forces  and  the  energy 
available  from  the  batteries.  The  test  bed  vehicle  was  used  to  validate  the 
synergistic  effect  of  all  models  in  EVSIM.  After  a  total  of  over  fifty  hours 
of  testing  the  EV  in  the  all  electric  and  hybird  models  in  both  urban  and 
cruise  speed  schedules,  the  EVSIM  program  was  found  to  be  within  ten  percent 


of  test  bed  empirical  data  in  all  parameter  values.  Additional  tests  were 
conducted  on  subsystems  such  as  batteries  in  static  tests  to  validate 
portions  of  submodels. 


AO 


V.  Additional  Concept  Investigation 


The  baseline  test  bed  vehicle  was  modified  to  investigate  various 
controller  schemes  and  control  laws  as  well  as  battery  charging  options  and 
vehicle  safety  operation. 


Control  Schemes 


While  the  energy  storage  devices  (batteries)  are  the  recognized 
foremost  deficiency  for  EV's,  the  power  conditioner /controller  seems  a  logical 
choice  to  investigate  improvements  in  both  performance  (drivability  and  com¬ 
patibility  with  existing  traffic  flow)  and  efficiency  (range). 

The  three  controllers  investigated  are  representative  of  the  types  of 
controllers  in  use  today  or  proposed  for  EV's.  Type  1  controller  is  an 
electro-mechanical  three  state  (discrete  energy  level)  contactor  controller. 
It  is  simple  and  reliable  in  its  operational  routine.  Figure  15  shows  a 
simplified  block  diagram  of  this  controller  with  a  microcomputer  inserted  in 
the  loop.  Type  2  controller  is  a  six  state  device  which  has  a  number  of  both 
power  electronics  and  electro-mechanical  components.  While  this  type  of  con¬ 
troller  obviously  adds  performance  flexibility,  it  is  at  the  cost  of  oper¬ 
ational  complexity.  The  type  3  controller  is  an  effectively  continuous 
(infinite  number  of  states)  controller  which  is  primarily  power  electronics 
in  its  component  configuration.  The  type  controller  is  implemented  with  256 
states  when  the  microcomputer  is  inserted  in  the  loop. 

All  of  these  controllers  were  first  individually  configured  in  a 
manual  mode  (driver  operated)  with  mechanical  linkages  and  discrete  elec¬ 
tronics  providing  actual  inputs  for  the  manual  driving  tests.  The  con¬ 
trollers  were  also  individually  fitted  with  interface  for  a  microcomputer  so 
that  it  could  be  installed  in  the  control  loop  of  each  type.  It  was  then 
possible  supply  inputs  to  the  controllers  from  the  processor  through 
appropriate  analog  interfaces.  For  the  automated  driving  tests,  the  vehicle 
operator  would  demand  a  velocity  and  the  processor  would  interpret  this  input 
and  convert  it  to  appropriate  controller  commands  for  each  type  controller. 
These  controllers  were  modeled  and  their  characteristics  incorporated  into 
EVSIM. 


Microcomputer 

The  availability  of  computationally  powerful  and  fast  single  board 
microprocessors  make  them  the  logical  device  to  implement  the  control  schemes 
A  Z80  based  system  was  selected  (Reference  17).  The  Z80  is  a  bus  oriented 
system  requiring  minimum  additional  control  logic  for  addressing  and  reset, 
as  indicated  in  Figure  16.  A  Z80  microcomputer  also  uses  a  single  phase  two 
megahertz  system  clock  to  synchronize  its  operation  and  a  single  five  volt  DC 
power  supply.  The  serial  input/output  controller  (SIO)  has  two  independently 
programmable  data  communication  channe Is .  These  channels  are  used  to  drive 
the  RS-232  interfaces  for  the  operators  terminal  and  the  software  development 
computer.  The  Z80  counter  timer  circuit  (CTC)  contains  four  independently 
programmable  channels,  two  of  which  are  used  as  baud  rate  clocks  for  the 
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serial  ports  of  the  SIO.  The  two  remaining  counter  timer  channels  were 
cascaded  to  produce  a  104  Hertz  low  frequency  clock  and  a  eight  Hertz  clock 
for  the  control  system  time  reference.  Sample  period  is  used  to  generate 
interrupt  service  requests  eight  times  a  second  and  is  used  by  the  control 
software.  The  microcomputer  uses  three  parallel  input/output  controllers 
(PIO)  giving  a  total  of  six  8-bit  parallel  data  ports  for  use  as  input  or 
output  interfaces.  Each  PIO  port  is  independently  programmable  and  is  TTL- 
compatible.  When  used  as  an  output  port,  the  PIO  holds  its  most  recent 
output  state  providing  the  hold  mechanism  for  control  interfaces.  The 
microcomputer  also  employs  an  analog  data  acquisition  system  as  an  integral 
part  of  its  design  to  enhance  system  flexibility.  The  analog  data  acqui¬ 
sition  system  provides  the  common  functions  associated  with  the  collection 
of  analog  measurements.  It  also  uses  a  one  Megahertz  clock  to  operate  its 
analog  to  digital  converter  and  is  designed  to  be  interfaced  with  most 
microprocessors.  The  microcomputer  uses  two  types  of  memory  EPROM  and  RAM. 
The  EPROM  is  used  to  store  the  system  monitor,  application  programs,  and 
provides  a  means  of  correcting  software  errors  quickly.  A  single  power 
supply  (+5V)  EPROM  was  selected  to  minimize  power  source  requirements.  RAM 
memory  is  used  for  scratch-pad  storage  and  calculations.  The  RAM  is  also 
used  during  the  software  development  process  to  hold  new  programs  for  testing 
and  debugging.  The  low  power  Schottky  TTL  components  provide  the  address 
decoding  for  the  computer  memory  and  input/output  devices.  The  also  buffer 
the  Z80  central  processing  unit  and  are  use  in  the  reset  logic  that  initial¬ 
izes  the  microcomputer.  The  control  system  design  has  digital  interfaces 
which  translate  TTL  logic  level  signals  into  voltage  and  current  signals  to 
direct  the  power  electronics  in  the  high  power  portion  of  the  controller. 

These  interfaces  also  provide  the  power  supply,  signal  ground,  and  chassis 
ground  isolation  for  the  analog  data  acquisition  system  by  using  optical 
isolators . 

The  microcomputer  was  installed  in  a  research  and  development  mode  as 
shown  in  Figure  17.  Software  development  was  accomplished  quickly  and 
efficiently  in  a  second  more  sophisticated  microcomputer.  Upon  verification, 
the  software  was  programmed  into  EPROM's  which  were  subsequently  installed 
into  the  microcomputer  controller.  The  software  flow  is  shown  in  Figure  18. 
The  control  system  had  two  software  programs  stored  in  EPROM's.  The  first 
program  was  the  standard  "ROM  Monitor"  operating  system  for  the  microcomputer. 
The  second  program  was  application  software  which  implemented  the  control  and 
energy  monitoring  functions  of  the  control  system.  The  control  functions 
regulated  the  vehicle  and  engine  speed.  The  monitoring  function  kept  a 
running  value  of  battery  energy  consumed  during  a  driving  cycle.  Because  of 
the  need  for  positive  regulation  of  the  control  system  by  the  operator  a 
"foreground  and  background"  program  execution  strategy  was  used.  The  monitor 
operating  system  cycled  in  the  "foreground"  waiting  for  either  an  operator 
input  or  a  sample  period  interrupt.  The  application  program  was  executed  in 
the  "background"  on  every  interrupt  cycle  to  update  control  and  monitoring 
tasks.  The  monitor  also  contained  application  program  commands  to  initiate 
the  driving  cycle  and  report  energy  consumption  to  the  operator. 


Controller  Implementation  and  Results 


To  test  the  various  types  of  controllers,  three  techniques  were  used 
as  follows:  (a)  a  digital  computer  simulation  (b)  a  vehicle  test  without 
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Fig.  18  Software  Flow  Diagram  for  the  Control  Systt 
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microcomputer  control  (c)  a  vehicle  test  with  microcomputer  control. 


The  controller  schemes  discussed  above  were  implemented  first  in 
EVSIM  to  test  the  feasibility  of  each  type.  Both  the  urban  driving  cycle 
FUDS  (Reference  10)  and  a  constant  speed  cruise  segment  (88  km/hr)  were 
simultated  for  each  type. 

The  test  bed  vehicle  was  modified  for  the  road  test.  The  vehicle  was 
configured  with  the  eight  bit  single  board  microcomputer  described  above. 

Two  sets  of  tests  were  run  for  all  the  controller  types.  In  the  first  set, 
all  the  controllers  were  manually  operator  controlled,  and  in  the  second  set 
of  tests,  the  controllers  were  managed  by  the  microcomputer.  In  the  manual 
tests,  the  microcomputer  was  used  only  to  perform  some  safety  functions  and 
the  energy  monitoring.  In  the  automated  tests,  the  microcomputer  was  used  in 
the  control  loop  to  provide  velocity  feedback  and  filter  the  system  control 
commands  by  a  set  of  acceleration  limits.  The  microcomputer  used  stored  data 
and  a  subsequent  table  look  up  to  implement  type  1  and  2  cQntrol  laws,  and 
solved  an  equation  to  implement  the  type  3  control  law.  The  vehicle  tests 
were  conducted  for  both  urban  and  cruise  segments. 

The  results  of  tnese  runs  are  summarized  in  Table  IV.  The  energy 
used  and  efficiency  data  was  converted  to  an  efficiency  measure,  percent  of 
range  (km)  change  which  is  a  significant  factor  for  EV's.  The  digital  simu¬ 
lations  indicated  that  a  digital  computer  was  compatible  with  all  controller 
types,  and  that  the  efficiency  measure  increased  as  expected  with  the  number 
of  control  states  per  controller  for  all 

TABLE  IV 

Simulation/Vehicle  Test  Results 


Efficiency  Measure  (%  Range  Change)  Speed  Schedule 
Urban  (FUDS)  Cruise  Variation  (FUDS) 

_ (j,  km/hr) 


A-D 

0 

0 

9 

A-E 

7 

-3 

6 

A-F 

10 

-2 

4 

B-D 

4 

2 

7 

B-E 

10 

1 

4 

B-F 

13 

2 

2 

C-D 

7 

6 

6 

C-E 

14 

7 

3 

C-F 

19 

8 

1 

A  =  Manual  Operator  Control 
B  =  N 1 ;rocomputer  Control 
C  =  Simulation 


D  =  3  State  Controller  (type  1) 

E  =  6  State  Controller  (type  2) 

F  =  256  State  Controller  (type  3) 
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simulaced  driving  regimes.  The  microcomputer  based  controllers  showed 
efficiency  measure  increases  of  up  to  five  percent  range  (two  km)  over  the 
corresponding  manual  controllers  for  the  urban  driving  tests.  The  data 
shows  that  as  the  number  of  control  states  gets  large,  there  is  a  rapidly 
decreasing  marginal  increase  in  range.  The  microcomputer  based  controllers 
all  showed  a  consistent  two  km/hr  reduction  over  the  manual  controllers  in 
speed  schedule  variation  which  translates  to  better  traffic  flow  blending 
for  the  EV.  All  the  urban  driving  tests  verified  the  expected  result  that 
the  range  increased  directly  with  the  number  of  controller  states.  The 
cruise  segment  vehicle  tests  show  that  the  kind  of  controller  has  more  effect 
on  range  than  the  number  of  states  for  the  controller.  The  maximum  increase 
in  range  due  to  the  microcomputer  controllers  in  the  cruise  driving  tests  was 
two  km.  The  speed  schedule  variations  for  the  cruise  segment  were  numeri¬ 
cally  insignificant.  Additionally,  results  indicated  that  the  speed  and 
precision  provided  by  the  microcomputer  would  increase  the  efficiency  measure 
for  all  controller  types  compared  to  manual  operation  for  all  driving 
regimes. 

These  results  show  that  currently  available  EV  electro-mechanical  and 
power  electronics  configured  controllers  can  be  successfully  managed  by  a 
single-board  microcomputer.  While  these  microprocessor  based  controllers 
increased  the  efficiency  (range)  of  the  EV  tested,  these  increases  were  not 
of  the  magnitude  expected.  Besides  managing  the  controllers,  the  microcom¬ 
puter's  presence  in  the  EV  was  found  to  be  very  convenient  and  useful  as  it 
is  in  a  growing  number  of  conventional  petroleum  fueled  vehicles.  The  micro¬ 
computer  performed  safety  functions  such  as  monitoring  the  EV  systems, 
checking  for  malfunctions,  and  served  as  a  practical  way  to  monitor  energy 
usage.  This  energy  monitoring  provided  an  accurate  state-of-charge  indicator 
for  the  EV.  The  microcomputer  could  then  act  as  a  manager  for  battery  charg¬ 
ing  as  well  as  normal  discharging  which  will  provide  for  improved  battery 
operational  performance  and  long  battery  life. 

Optimal  Control  Law 

Because  of  the  improvement  made  by  the  multi-state  controller  (type  3) 
just  discussed,  it  was  logical  to  investigate  whether  an  optimal  control  law 
applied  to  this  plant  would  produce  significant  results  in  terms  of  both 
efficiency  (minimize  energy  used)  and  performance  (minimize  deviation  from 
demanded  speed).  The  first  improvement  results  directly  in  more  range,  and 
the  latter  improvement  results  in  increased  acceleration  performance  or 
drivability . 

Optimal  Control  Design 

The  optimal  control  based  controller  developed  in  this  effort  was  of 
low  order  to  match  the  reduced  order  model  of  the  EV  used.  The  reduced  order 
model  shown  in  Figure  19  was  appropriate  for  this  effort  since  its  accuracy 
has  been  shown  to  be  high,  and  the  purpose  of  the  optimal  control  law  investi¬ 
gation  was  to  see  if  the  magnitude  of  the  potential  for  improved  efficiency 
and  performance  justifies  further  pursuit. 
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The  dynamics  of  the  EV  system  shown  in  Figure  !9  allow  several  further 
simplifications.  The  speed  of  the  digital  processor  and  associated  elec¬ 
tronics  make  their  effects  negligible  when  compared  to  the  drive  train  com¬ 
ponent  responses.  Similarly,  the  time  delay  (t^)  is  negligible  compared  to 
the  dynamics  of  the  other  components  in  the  power  controller.  Since  an 
automated  controller  is  sought,  a  set  point  controller  seems  a  logical  choice. 
This  controller  would  tend  to  minimize  the  effects  of  disturbances  (wind, 
road  grade  and  roughness,  etc.)  and  still  allow  an  operator  interaction  to 
vary  the  set  point  for  speed  schedule  changes  and  safety  considerations.  The 
resulting  simplified  set  point  controller  is  shown  in  Figure  20.  The  result¬ 
ing  reduced  order  state  space  model  is  as  follows: 


(60) 


(61) 


and 

Y  =  CTx  =  [o  l]  x  =  Va 


x  =  Ax  +  bu+£d 


where  x^  is  the  motor  voltage,  V  ,  x^  is  the  vehicle  velocity,  V  ,  d  is  the 
disturbance,  and  the  input  u  is  ?*he  controller  voltage,  V  .  Theatime  response 
of  this  system  with  no  feedback  for  a  unit  step  input  forCV  and  zero  dis¬ 
turbance  is  shown  in  Figure  21.  As  can  be  seen,  the  systemChas  a  stable 
overdamped  response  which  illuminates  the  performance  shortfalls  of  this  type 
of  vehicle. 

For  this  set  point  controller,  the  quadradic  performance  index  is  as 

follows : 

J  =  4/(x-x  )TQ(x-x  )  +  r  (V  -V  )2dt  (63) 

- s - S  CCS 

where  V  is  the  set  point  controller  voltage  required  by  the  system  to  reach 
x„  in  steady  state,  =  £xjg,  1  ;  x.  is  the  motor  voltage  at  set  point 
velocity;  x^s  is  the  set  point  velocity  letermined  by  the  operator  or  speed 
schedule,  J  is  the  cost  functional  to  be  minimized;  Q  is  the  weighting  matrix 
penalizing  deviation  from  set  point  state  values;  and  r  is  the  weighting 
scalar  penalizing  the  use  of  control.  For  this  linear  quadradic  set  point 
regulator  formulation,  an  optimal  control  V  *  can  be  found  which  will  mini¬ 
mize  the  cost  functional  J.  Since  both  stages  x  and  x»  are  controllable  by 
the  input  V^,  and  both  are  observable,  the  problem  is  properly  posed  for  full 
state  feedback  optimal  control  (Reference  12).  The  optimal  control  law 
which  minimizes  the  cost  functional  of  Equation  63  is: 


V  * 
c 


+  V 


cs 


(64) 
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where 


K*  =  (  1/r)  b‘P 


The  matrix  P  is  the  steady  state  Ricatti  equation  solution  for  the  linear 
quadradic  set  point  problem. 

ATP  +  PA  +  Q  -P  b  bT  P  (j)  =  0  < 

The  set  point  is  determined  through  the  scheduled  vehicle  velocity,  V 

3  S 

That  is 


- — V 

K2  as 


0.8  V 


1.0  V 


V  =  *  U—  V  =  0.13  V 

cs  Kj  K2  as  as 


The  Q  weighting  matrix  which  penalizes  state  deviation  from  the  set 
point  can  be  used  to  tune  the  optimal  controller  toward  higher  performance 
(quicker  response,  better  drivability ) .  The  r  scalar  (using  a  single  control 
input)  which  penalizes  the  use  of  control  can  be  used  to  tune  the  optimal 
controller  for  increased  efficiency.  As  expected  with  any  quadradic  per¬ 
formance  index,  the  two  objectives  of  increased  performance  (through  large 
values  in  Q)  and  increased  efficiency  (through  a  large  value  of  r)  drive  the 
optimal  feedback  model  in  opposite  directions.  Because  there  is  no  single 
combination  of  Q  and  r  that  produces  optimal  performance  for  both  objectives, 
several  compromise  combinations  of  the  weighting  elements  were  investigated. 
The  resulting  closed  loop  system  with  the  optimal  feedback  gains  is  shown  in 
Figure  22.  The  gain  K  *  is  used  to  adjust  the  final  output  for  correct 
steady  state  response.  Only  diagonal  elements  of  the  Q  matrix  were  varied 


TABLE  V 


Weighting  Parameter  Variations 


Case  1 

1 

1 

1 

Case  2 

1 

10 

1 

Case  3 

1 

100 

1 

Case  4 

1 

1000 

1 

Case  5 

! 

1 

10 

Case  6 

1 

1 

100 

Case  7 

1 

1 

1000 

TABLE  VI 


Closed  Loop  Eigenvalues  and  Feedback.  Gains 


Eigenvalues 

V 

K  * 

2 

V 

Case 

1 

-60.82 

-0.79 

0.852 

0.475 

9.70 

Case 

2 

-60.79 

-0.20 

0.872 

2.438 

24 .70 

Case 

3 

-60.50 

-6.26 

0.937 

9.155 

75.85 

Case 

4 

-57.11 

-20.93 

1.125 

30.607 

239.  14 

Case 

5 

-21.44 

-0.75 

0.  194 

0.  137 

3.22 

Case 

6 

-11.66 

-0.60 

0.029 

0.028 

1.39 

Case 

7 

-10.  18 

-0.51 

0.003 

0.003 

1.04 

Case 

8 

-  5.25 

±j3.86 

1 .0 

0.0 

8.59 

K^*  assures  appropriate  steady  state  response 
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)  because  it  was  evident  the  q0„  element  corresponding  to 


state  x„  (velocity)  had  the  main  influence  on  improvements  in  drivability. 
Controllers  were  designed  for  the  set  of  Q  and  r  elements  shown  in  Table  V. 

Cases  2  through  4  represent  systems  optimized  for  performance  and 
Cases  1  and  5  through  7  represent  systems  optimized  for  efficiency.  No 
cases  of  either  r  =  0,  or  both  and  q  *  0  are  included  because  only 
physically  meaningful  and  realizable  systems  were  sought  (Reference  24).  The 
effective  closed  loop  eigenvalues  for  the  cases  are  shown  in  Table  II  along 
with  the  feedback  gains.  Case  1  for  the  efficiency  oriented  systems  seems  to 
be  a  good  compromise  between  speed  of  response  and  efficiency  achieved.  Case 
1  was  chosen  as  the  efficiency  candidate  to  investigate  further.  Case  3  was 
chosen  for  further  investigation  from  among  the  performance  oriented  systems 
because  of  its  compromise  between  speed  of  response  and  performance  in  terms 
of  minimized  state  deviation.  Case  8  (Table  VI)  simulates  an  operator  in  the 
feedback  loop  modeled  as  a  pure  gain  (k  =  1). 


Optimal  Control  System  Implementation  and  Testin 


The  systems  represented  by  Cases  1,  3,  and  8  were  implemented  in  the 
electric  vehicle  digital  computer  simulation  EVSIM.  The  controller  sub¬ 
routine  in  EVSIM  was  modified  to  incorporate  the  optimal  control  laws  of 
Cases  1  and  3  as  well  as  the  control  law  for  Case  8.  Figure  23  shows  the 
comparative  system  responses  to  a  step  change  in  set  point. 

A  separate  look  at  the  various  systems  responses  to d isturbances  shows 
Cases  3  and  8  respond  about  equally  well  to  compensate  for  road  grade  and 
wind  disturbances  for  satisfactory  performance.  The  response  of  Case  1  to 
disturbances  demonstrates  a  lack  of  acceleration  performance.  Figure  24 
shows  the  relative  system  responses  to  a  simulated  hill. 


System  Response  Comparison 


The  results  of  the  simulations  are  displayed  in  Table  VII.  Comparing 
Cas^s  1  and  3  in  the  driving  segment,  it  can  be  seen  that  Case  1  produces  a 
much  more  efficient  system  but  at  considerable  reduction  in  acceleration  per¬ 
formance.  Case  3  dominates  Case  1  in  performance  but  at  a  significant  cost 
in  efficiency.  The  relative  separation  between  the  two  cases  is  less  for  the 
cruise  segment  but  the  trends  are  preserved.  The  overall  performance  of  Case 
8  was  surprising  in  that  it  came  close  to  matching  the  drivability  of  Case  3 
and  the  efficiency  of  Case  1  for  both  the  driving  and  cruise  segments.  Com¬ 
pared  to  manual  controllers,  Cases  1  and  8  produced  approximately  33  percent 
increase  in  range  capability  while  Case  3  produced  a  24  percent  decrease. 

The  corresponding  speed  schedule  variation  was  improved  by  approximately  132 
percent  for  Cases  3  and  8,  and  decreased  by  12  percent  by  Case  1.  Compared 
to  conventional  microprocess  driven  automated  controllers.  Cases  1  and  8 
increased  range  capability  by  approximately  11  percent  while  Case  3  showed  a 
decrease  of  35  percent.  The  corresponding  speed  schedule  variation  was 
improved  by  approximately  13  percent  for  Cases  3  and  8  and  decreased  by  34 
percent  for  Case  1. 

TABLE  VII 

Control  System  Comparison 
Five  Minute  Driving  Segment 
(No  disturbances) 

Energy  Used  RMS  Speed  Schedule 


(kwh ) 

Variation  (km/h) 

Case 

1 

(Efficiency) 

0.91 

3.84 

Case 

3 

(Performance ) 

4.35 

1 .  17 

Case 

8 

(Operator  in 
Loop) 

1.06 

1.90 

Five  Minute 

Cruise  Segment 

(With  disturbances) 

Case 

1 

1.20 

2.64 

Case 

3 

3.76 

.96 

Case 

8 

1.76 

1.28 

The  results  show  that  optimal  control  does  have  a  significant  effect 
on  EV  performance  when  compared  to  conventional  manual  controllers  but  that 
no  one  law  was  able  to  produce  improvements  in  both  efficiency  (range)  and 
performance  (drivability).  Implementation  of  the  optimal  control  laws  tuned 
for  efficiency  (Case  1)  and  acceleration  performance  (Case  3)  were  able  to 
produce  II  and  13  percent  improvement  respectively  in  range  and  speed  schedule 
variation  over  conventional  microprocessor  driven  automated  controllers. 

While  these  gains  are  notable,  they  are  not  of  the  significant  magnitude 
expected . 


The  performance  of  microprocessor  driven  automatic  controller  with 
an  ideal  operator  in  the  feedback  loop  (Case  8)  was  surprising  that  it  per¬ 
formed  well  against  all  optimal  control  laws  investigated  with  respect  to 
both  efficiency  and  performance.  This  would  indicate  that  simple  con¬ 
trollers  may  be  able  to  produce  satisfactory  results  for  an  EV  class  of 
systems,  and  that  as  this  effort  indicates,  more  complex  controllers  add 
only  marginal  improvements  to  system  operation. 

Cases  3  and  8  performed  satisfactorily  in  handling  disturbances 
such  as  wind  and  road  grade  changes,  but  neither  was  tuned  for  this  purpose 
It  is  expected  that  a  system  specifically  tuned  to  minimize  the  adverse 
effects  of  these  disturbances  would  increase  drivability  and  thereby 
increase  the  attractiveness  of  EV's.  Such  an  area  seems  potentially 
attractive  for  future  efforts. 


Summary  and  Conclusions 


The  electric  vehicle  simulation  program  EVSIM  does  satisfy  the 
primary  objective  of  providing  a  useful  evaluation  and  design  tool  for 
electric  vehicle  systems.  The  program  has  the  flexibility  to  simulate  both 
electric  and  electric/hybrid  configurations  using  series  or  shunt  motors, 
step-voltage  or  continuous  controllers,  and  gear-ratio  or  continuously  vari¬ 
able  transmissions  for  any  general  shape  and  size  of  vehicle.  Any  specific 
driving  cycle  may  be  used  to  simulate  the  vehicle  operation.  EVSIM  has  the 
ability  to  internally  modify  a  basic  input  vehicle  to  provide  improved  per¬ 
formance  to  meet  required  standards. 

The  general  conclusion  formed  during  the  validation  of  EVSIM  was 
that  it  did  indeed  provide  an  easily  obtained,  accurate  evaluation  of  a 
vehicle's  potential  performance.  With  very  little  computer  time,  it  clearly 
illuminated  the  large  and  small  effects  of  vehicle  component  modification. 
This  ability  should  certainly  make  the  program  an  efficient  tool  for  both 
the  designer  and  vehicle  evaluator. 

Additionally,  hybrid  power  augmentation  is  required  to  meet  the  per¬ 
formance  standards  set  in  Chapter  II  for  general  purpose  USAF  passenger  cars 
The  energy  available  from  the  state-of-the-art  lead-acid  batteries  is  simply 
insufficient  to  provide  reasonable  high-speed  cruising  range  when  con¬ 
strained  to  the  dimensions  of  a  compact  vehicle.  Until  significantly 
improved  batteries  become  available,  a  parallel  hybrid/electric  vehicle  is 
the  most  feasible  alternative  to  the  gas-powered  conventional  USAF  sedan. 

It  is  not  clear  that  the  increased  efficiency  afforded  by  a  micro¬ 
computer  based  controller  justifies  the  microcomputer  for  EV’s.  But  this 
increased  efficiency  capability  coupled  with  a  decrease  in  traffic  flow 
speed  variance,  and  the  availability  of  an  accurate  energy  usage  and  safety 
monitor  suggest  that  the  microcomputer  will  make  the  EV  technically  and 
economically  feasible  for  the  average  urban  commuter. 

While  increases  in  efficiency  and  performance  can  be  gained  by  the 
application  of  an  optimal  control  law,  these  gains  were  not  significant  to 
overcome  deficiencies  presented  by  the  batteries. 
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EVSIM  program  listing 


PRC  C  9 1  ■  E  V  $  I  »  (  IaP‘jT,!3lTP>jT.Tap£$-|NPUT.Ta;>E6-CUTPUT> 

Cl  »£AS ION  YX(l<1CG),CB<ICIiCSTP<SI,»STPIlQI»ETA*MlO«3C),ETATIlQ«2QI 
CtT^TPtl0».'J8AU5l,»SiU|Sl,V«lv.0O) 

REAL  wi«,w,>,S.»<>,CT,NSTP,-l,',CF,HSH 
FC«"AT(  tCF«.l  t 
fcrpati  i8f  a.  1 > 

P  P  I  A  T  I  ft  ,  <3  C  2  I 

F cap* 1 < in  1 , 3  1  X 7H*  •«EL£CT*IC  VEHICLE  S  l  *  U 

CLATIOA«**| 


EATER  VEHICLE  CONSTANTS. 
HASStPwl,  FIxEC  "ASSCFwt 


CRAG  CQEFF (CO  I  .FRONTAL  A  8  E  A ( F  A  )  t  PAYLQAQ 


»EAC».CO,FA ,PP  ,FH 
PR IATIA ,QC3  I 

FC8PAt(///2Cx22HVEHlCL6  CONSTANTS  ARE:! 

PfltNTtft,9CAICC.FA,PW,FH 

FCRPAT<35»t7HC»AC  COEFFICIENT-, F5. 2, /35X13HFRCNTAL  AREA-, F5. 2.1  AH 
CSCUARE  w£ TER S ,/ 35X1 Thpaylqac  P A S S • , F A .0 , 1  CM  K I LOC P A P S  , / 35X UHF  I  XE 0 
C  PASS-, FA. 0» 

EATER  ENVIRONMENTAL  Constants.  apbienT  T E p PE 8 A TU8 E I T PP I »  ELECTRIC 
POWER  COS  T ( ECC  ST  1 ,  PETROLEu-  FUEL  COSTIFCCST). 

REAC*.TPP,ECCST,FCOST 
PR  I  A  TI5.905  » 

FCRPAT  (  /  70X28fEavircnpEATal  CONSTANTS  ARE:) 
PRIATIfc.RCfclTPP.ECOST.FCQST 

FCR"ATI3SX2CHAPRIENT  TEP?ERAT(.RE-.FA,0,aH  CELSIUS. /35X23HELECTRIC 
C  £  A  £  R  G  V  COST-  1,Fa.3,9h  P£8  w -h H , / 3 5 X 2 2  HP E T R CL E UN  FUEL  COST-  S.FA.2 

c. ich  per  Linn 

Enter  VEHICLE  power  REDUCTION  VARIA8LES.  ROLLINC  8A0IUS1RR),  TIRE 
ROLLING  COEFFKIENTICRI,  a  t  L  E  8AT10(aR»,  nu.-BER  OF  FCRwARO  SPEEDS 
<0  IF  CVOUSl.  .NAX  VEHICLE  CESIGN  SPEEOIVNAXI. 

RE*C*.8R , C  R , A  S , I S , VP A  X 
PRINT(6,9C7> 

FC8PATI/20X30EP0WER  REOUCTICN  VARIABLES  ARE:) 
PRIATC6,RCB)RR,CR.A8,tS,VPAX 

FC***T(3S*2CH1  [RE  ROLLING  R a C l US- ,Ffc. 1 , 7H  PE T E R S , / 3 5 X 20mR OLL I NG  CO 
CEFFICIEM-.F  7.A./3SXIU  a/LE  R a T 1 P • . F 5 . 2 » / 3 5 X 2 1 HNU " B E 8  F  0  RWA8  0  GEAR 
CS- t  12 ,/3S» I ?HP A  X  I *UW  VELOC I T Y-,F A.0.6H  AH/HR) 

IFITS.NE. 0)10.20 
CONTINUE 

EATER  TRANSPISSION  RATICSXIF  »ULT I S PE EO > I CR I  I S  I )  ON  FIRST  CARO.  ON 
NEXT  CARO  EATER  OESIREO  NOTCH  UPSHIFT  AnQ  DOWNSHIFT  8 P N ( UP  SHF T , 

ON  SHF  T  )  . 


REAC*,ICRI|),I-I,1S) 

PEAC«,UPSHFT,CASHFT 
PR  I  A  T (6  «9C9  ) 

FCRPAT(/20X3<THF  TXEO  SPEED  TRANS"1SSI0N  VARIABLES  ARE:) 
PR|ATI6.9)0)LPSHFT,CNShFT,C8ll) 

FCRPATI 3SX2CUPCTCR  UPSHIFT  S P E E 0 - , F S . 0 , AH  8 P H  ,  /  35X 22 H*0 T OR  OOVNSMI 
CFT  SPEEC-,FS.C,ah  RP«,/3SX12HCEAR  R A T 1 0 S • , F6 . J  ) 


if  (is.Eo.mc.ii 
11  00  I  2  1-2,  is 

60  12  PR  I  n  » 1 1  ,911 1 C a  (  '  ) 

811  fC»HII  W»,f6.3l 
CC  10  30 
20  CONTINUE 
C 

t5  C  ENTER  C  V  T  VAR|JRLES.  RATEO  INPUT  POkERTPRATTI  ANO  Nl-BFR  CF  SPEED 

C  RATIO  STEPS  £  C  p  INPUT  CA  T  A  C  s  c  I  3  S  1  CARO,  E  E  E  I  C  I  F  n  C  Y  QAlA(ETAT) 

c  ON  M»r  C  A  8  C  SfSIES:  F  I  3 '  r  CAPO-  S  P  £  ?  0  PATIO  STEPSTTSTP),  FOL- 

C  t-  C  i  C  C  A  8  C  S-  £TAT'E.CT.O  A  NC  .LT.1>  IN  10F5.1  FOR  EACH  SZ  OF  PRATT 

C  FROM  SZ  T01CCZ.  2  C  A  8  C  S  8£0C  FOR  EACH  TSTP>. 

7C  C 

RE  A0»,  PRATT  ,  ITS  TP 
R£AC».(TSTP(U,I-l,ITSrP) 

OC  21  I-1.I1S1P 
.  00  21  2-1.2 

7«  21  RE*C(S.8C0ME1AT«I.(J-H*10«k*,k-1,10I 

PRTM(t,9l2» 

912  FCRmAI(/20XA9NCCNTINUCUSLY  VARIABLE  TRANSMISSION  VARIABLES  ARE:) 

PR  IM<6.913  )PSA  TT 

913  FCRPAT(3SXl3HPCh£R  8 A T I NC - , F 5 . 1 , 3H  KW) 

8C  30  CONTINUE 

c 

c  ENTER  MOTOR  VARIABLES.  "0T03  T Y P E I » T r P » < 0- Shun T , l • S E R  I E S > .  K  FAC- 

C  TC»|KPOT),  A  8  *  A  TUP  E  R £ S I S T A N C£ ( 8  A R N I .  FIELD  R £ S I S T A N C E I R FL 0  I , 

c  RATEC  POwFR I  PS i THJ  ,  max  Prw  -LL  T  (  I  P  a  ax  )  ,  max  CuRTA-ax),  BASE  SPEEO 

es  c  cespoi.  anc  m^ber  cf  speec  steps  fc«  in.>ut  oata  on  first  caro. 

C  ENTER  £PF  IC  IEN  C  v  f  E  T  A-)  CN  *. e  x  T  CARO  SERIES:  FIRST  CARO-  MOTOR 

C  SPEEO  RATIO  STEPS<8pM)(*Sr?l  •  FOLLOWING  CAPOS-  ETAM<.CT.O  AND 

C  • L  T . I >  IN  ICFi.l  FOR  EACH  ICZ  OF  PRATN  UP  TO  30QZ<3  CAROS  8E00 

C  FCR  EACH  MSTPJ. 

SC  c 

REAC«,MTYPt*MCT,RARM,RFL0,PRATM,IPMAX,AMAX,8SP0,lPSTP 
R£AC*,(MSTP( 11,1-1, IMSTP) 

OC  31  I-I.IPSTP 
00  31  J-t.3 

9«  31  REACTS.  SOOHETAMtl, »  J- 1 1  ♦  10  ♦  K  I ,  x  .  i ,  10  » 

PRTNTT6.9MI 

9  1  A  FCRMATT /2CX20FMOTQR  VARIABLES  ARES) 

IFTPTTP. EC. 0132,33 

32  CCN1INUE 

IOC  PRINTI6.91SI 

915  FCRPATT35X17HPCTOR  TYPE-  SHUNT) 

GO  TO  3A 

33  PRINTT6.916) 

916  F0RMATI3SX1EHPCTCR  TYPE-  SERIES) 

1C5  3A  PR  IN T Tt ,91 7  I  PS  a f m.BSPO , A-ax 

917  FCRPATUSX12KATE0  POv£R-,F5.1,3H  K  w  ,  /  3  5  X  l  1H  B  A  S  E  S  P  E  EO  •  ,  F5 . 0  ,  AH  RP 
CM./3SX21HMAX  ARMATURE  CURREN r. ,FS.0,5H  AMPS) 

AO  CONTINUE 

c 

110  C  ENTER  CONTROLLER  VARIABLES.  ENTER  TYPETCTYP)  CN  FIRST  CARO:  l  IF 

C  FINITE  STEP  SNITCHING,  2  IF  CHOPPER. 

c 

REACA.CTYP 
IFTCTTP. EC. 1.1*1,50 


CONTINUE 


1  1 ?  <.1 
C 

c 

c 

c 

120  C 


12*  6  2 
920 

9  21 

1  3  C  90 
C 

c 

C' 

c 

13? 

929 


1H 

60 

61 

C 

C 

16?  C 

c 


15C  926 

C 

c 

c 

15?  C 
C 


I6C 


927 

16? 


70 

71 

7C  C 
C 


EVIfU  FINITE  STEP  CCNTCCLLf*  VARIABLES.  F  !  R  5  T  CARD-  NL-BER  OF  VCLT 
ACE  ST  E  P  r.  I  I  CS  1  ?  I  .  NEXT  CARO-  fBACTICN  Of  I  AL  VOLTAGE  ECU  EACH 
STEPICST?!,  A  a  C  EffECTIVE  SCuiCE  R  £  S  I  S  T  a  nc  E  <  B  9  a  T  »  F0»  EACH  STEP. 

B  £  A c  *  »  ICSTP 

»£Ar«,fcsTPTn.t-i,icsTpi,(BNiT||i,i.ltlcsTpi 
PS |NT (6,920)  ICSTP, C$T?(  1  |  «B?AT(  1  ) 

00  62  1-2. ICSTP 

PC  I N  T  I  6 ,92  l  )CS  IP  tl)  . »9 A  T  I  1  1 

ECS*aTI/2C*I2.21h  STEP  CONTROLLER  VARIABLES  A  8 E : . / 3 5  * l 7H VOLT  AC E  FR 
CAC T ICA • ,f S . 3  ,  3 V  I  7HS OLSCE  l -P E 0 anCE • .F a . a , 5M  OH-ST 
F0R-ATI52XF5.3.20vF5. 31 
CO  TO  60 
CONTINUE 

ENTER  CHOPPER  CONTROLLER  SOLACE  I -PE 0 AnCE  I  R B A Tl  .  MODEL  ASSUMES  1.5 
VCLTS  JUNCTION  LOSS  ANQ  37  RESISTIVE  LOSS. 

REAC*.R0AT( 1 1 
PR  INTI  6 .929TRB ATCl » 

ECS»ATI/2CX33fCHCPPER  C  CN  T  8  CLLE  8  VARU9LES  AR E : » / 3 5 X l 7HS OURCE  tMPE 
CCANCE- .F5.3 .5H  CH-S./35X.32HASSU-E0  JUNCTION  CROP-  1.5  VOLTS, /35X2 
CfTHASSU-EO  RESISTIVE  LOSSES*  3  21 
CSTPIll-l.  1  ICSTP-l 
OC  tl  1-1  .ICSTP 
R3ATII  tl-RBAK  II 

ENTER  PATTEST  VARIABLES.  FIRST  CARO-  NO-INAL  TOTAL  V OL T A CE I VNOM  I 
ANC  TOTAL  MASSIB-I. 

REAOA.VNOM.BM 
0  P  1  •  9  P 

PR  INTT6.9261VNC-.PM 

FCR*ATT/20*22F?attE*V  VARIABLES  ARE :  , / 3 5* . 1 6HN0- I N aL  VOL T ACE • . F 6 .0 
C»/35*13H BATTERY  -ASS-.F5.0I 

ENTER  HVBRIO  VARIABLES.  FIRST  CARO-  l  JF  HTBRIC.O  IF  NOT.  NEXT  CARO 
< I F  IHre-lJ-  FV8PI0  SPEEO  RATIOIHRI,  COUPLING  E F F I C  I  E N C T ( E T A H  )  , 
ENCACE-ENT  SPEEOIvhT,  ANO  NC«pal  OPERATING  TRANSMISSION  CEARIlSH). 

REACA.iHra 

HP-O. 

VH- IOC. 

IF! IHT8.EC.0IC0  TO  71 
REAC* .HR ,ETAH  ,VH, ISH 
MP-10. 

P»INTI6,927IHR,ETah,VH,ISH 

FCRPAT I/2CX  .  3  2  F  7*  w  SEBIES  HVpS|o  VARIABLES  AR E : . / 35 X 25HH V8 R I  0 / -0  TO 
C»  SPEEO  RATIC-.F6.3./35X2CHCCUPL1NC  E F t  I C I E SC T • , F 5 . 1 , 1H - , / 3 5 x  1  7he N 
CCACEPENT  SPEEC-.F5.1.5HXM/HR./  35X22HTRANSPISSICN  PO S l T I  ON- , 1 2 . 5H  C 
CEART 
IS-  ISH 
CONTINUE 

PROCRAM  IS  NOV  REA01NC  INPUTTED  ORIVlNC  CYCLE  SPEEO  SCHEDULE  OATA. 


C  F  !  a  *5  T  C*BC-  T  I  »  E  INO?  «F  *.  T  t  T  >  ,  Su*?EB  CF  OaTA  p  C  I  M  $  (  h  D  A  T  A  )  .  FOL- 

C  ICuHC  CAPOS-  SPEED  Ufi  I  i.L  E  i  v  t  |  pi  l-F-.l.  LAST  CABO-  COVERSICN 

C  FaCTC»<0  IF  vx  In  hPn,  1  IF  vx  In  a  bm  >  ( CCNy  | , 

C 

RE  AC*  ,T .NCATA 
NNC  » T  A -AC  A T A / 1 8 ♦ l 
DO  (Q  J-l.NNCATA 

80  BEAC«5.<;ClMVM18«IJ-ll*K»,K»i,l9| 

BE  aC  *  ,CCNV 
IF  I  CON V.EG.C.  I  El, 82 
ei  oc  e i c i  j*i»ncata 
«ici  v*i j»n«v* u ♦  n« l.602 

82  NANC A T A- ND AT  A / t C ♦ l 

EATER  PROGRAM  OPTION.  FIRST  CARO-  0  IF  SVSTE*  EVALUATION.  I  IF 
SvSTEP  C6SICN.  A£  x  T  CAROOF  IPRQC-1>-  NAXINUN  ALL  Cw  A8L6  BATTERY 
DISCHARGE COS-AX  I  ,  ano  ilXIACF  ALLOWABLE  SPEED  SCHEDULE  0EVJATI0N 

I0VPA1I. 


READ*.  IPHOC 

1FIIPR05.EC.CI8201.8202 
62C  1  PRlAT(6,<nOI 

<»3C  F0«»AT(/2OX33HSVSTE«  EVALUATION  OPTION  SEIECTEOI 
CC  TO  8203 

8202  HEAC«,0SP4X.0VPAX 

PRtATl6.O31ICSPAX.0VN AX 

<*31  FGRPAU/2CX3SFSYSTEP  OESICN  OPTION  SELECTEO  U  I  TH  :  /  35  *2  2H  N  A  x  BATTER 
CV  01SCHARC£*.F«,l,2M  JI./3SX2CHNAX  SPEEO  OE  V  I  A  T  I  ON  •  ,  F  8 . 1 , 6M  KN/MRI 
8 2C 3  OEV-O. 

OIS-O. 

83  VVl-O. 

vm-c. 
vx  i  n-o. 

U  •  0  . 

N8-C. 

NH«0. 

0-0. 

KK  ■  t 

PPAX.PRATP* IPPAX 
np*S.7S**4LCGIC(PRaTN|*21.5‘HN 
n»( pp«8"*FN 1/ .78*.83*PH 
AH  T  «0 . 

IFICIS.EO.O. 18 302,8301 
8301  FA-  (P/Pt ) «• .S*F A 
8 3C2  Nl-P  t  ClS*0. 

00  160  TI-l.AAAOATA 
PRIATI6.928) 

<128  FCRPATI  1H1 ,  13CH  TI«E  VELC-SCH  VELC-aCT  OEVIATTON  F-ROLL  F-AIR 
C  F-ACCEL  INPUT-POWER  gear  SHAFT-RPN  CONTROLLER  nTR  aRn 
C8AT  OISCHARCEI 

PRINT*,"  ISECI  IKN/MRI  (*-/f*R|  (K«/HR)  (Nil  (NT)  (NT 

C)  IPTRI  (Hvei  SELECT  INTRHhVSI  STEP/FIELO  VOLTS  APPS  (X 

Cl" 


PRINT*,"  • 

oc  itc  irr-i.ac 
I  •  I  I  I  ♦ I  I *60-88 
IFIVXHI.CE. <78. <71170, 85 


63 


nc 


231 


2*C 


2**. 


2 1C 


21! 


.  2fcC 


2*1 


27  C 


27! 


2  8C 


281 


81  V9AT-(73.l«.C2«9»3t1-.CClCl«C!1«0m»V*0‘*/72. 

U  *  *0 

90  V  \r  -  V  C  !  t  *  .  2  7  7  8 
CELIAV-.02 
A»0. 
vi  r.Q. 

*i  -o. 

P1H.0. 

ppp.Q 

P9CCH«0. 

PAVL«0. 

PCFM.02 

VVA-(VVl*VV|/2. 

(FM1.EC.0712.il 

*71  CALL  PTRSPOMVA.RR.AR.GB.UPIHFT.ONSMFT.IS.KK.KII 
«CF-*CF*.C0C2<MM.i9»GR<KK»/*8|#*2. 

12  FR*(CR*.CCCCCC2C9»vVa«»2.8)*p*9.81 

F8I«*VVA**?.«1?6.<.*FA*CC/I273.*T.*P) 

F  ACC»PM  W-V  V  I  l/T«.1CF 
F»FR  »F  A  JR  «F ICC 
P»F«VVA/1000. 

IFIVVA .EC.O.KC  TO  111 
9C02  P$P0«V| I MVAJMOO. 

PP.P/PNAXMCO. 

IF(PP.EC.0.(92C2.920l 

920  ETAA»95.8-I.01217*P1PCPP!P0«.8879PPSP0**.2611/ABSCPPI 
92  C  2  IF  Ml. EO. 0)920, 920* 

92  C 3  *1*eiP0 

CPU  TRMEFF(P,ETAA,?RATT.V*A,RR,AR,91P0.TlTP,tT5TP,ETAT,ETAI 
CC  TO  ICO 

9  20*  !F(PP.  EC.  0.  Mil.  9201 
92C1  P1PCT.P1PC*C»|PK  J/C8MS1 
IFICMKKI.LT. 1.193. 9* 

S3  E*A>>E*P<  .03*M.-l./CR(KKM  f 
CO  TO  99 

9*  E’AI’EXPt .03*11. >C8(KK|  M 

99  ETA-6PPXPI99.-I.  CO  <.C01*P1POT«P1POT*.19*9*P$  POT- 1.16 11/481  (PPM 

100  IFIPP.CT.C* 71C1.1C2 

101  PRCO*P/ETAA/ETA*lOOOO. 

co  to  Mic.uci.uancTvp 

102  P4VL*P*f T44«ET4/10000. 

CO  TO  111 

UO  C*lL  7*!TEP(P«C0,.H,B49P,K»QT,V9AT,tC1TP,ClTP,*,P*Vl,4P4X,PBAT*t 

CPSTP,(»lTP,ET4P,P»'*,»84T,VLT.tMlrd,HB.ETAH,P»QCP.»(I)*VM,KTTP,8FL01 
CC  TO  III 

1101  CALL  VCFCPIFTTP, BAB*. aFlO,KPGT,A«AX, PRAT*. I«STP,*lTP,ETA*,fi BAT, VBA 
CT.IFTe,NR,£TAFtVH,PKCG,',1fVM>,A,VLT.PAVl,PRC0H» 

111  lF(APSMP4VL-fBC0t/f/VMM.CT..000l)0ELTAV-.5 

112  IF  (kkk.EO.2M2  1 .118 

118  IF(PAVt.Cr.PRCCMl9,l20 

119  VII  J«VM»*CELTA» 

KK  K  •  l 

CO  TO  90 

120  IF(pkk. EC. 1)112, 121 

121  IF(PAVL.LT.PBCC)  122,112 

122  Vl|  )«VM  l-CELTAV 


KK*  -2 

CO  TO  <50 

HI  W?  •»  t  ♦  T/  3600  .  *P  AVL 
CO  TO  153 

152  w-a«T/16C0.*PAVL 
NN-AH»7/<6CC.*PR00H 

153  IT-M1-T 
VV  t  -VV 

PflCCF-PAVL-PRCCH 

IF<v<I).CE.VH.AN0.lWY8.E0.1l-5H.-S*HR 

o«o»w4»r/icoc. 
ne  v.vm- vx  in 

PR|MC6,9iaitT,Vxm.V<n,CEV.FR.FAtR.FACC.PRaOE,PRCOH,KK.NS.R5H,* 

C--.FLC.VLT, A, C  15 

918  FCfl"ATIlXt5»FS.l,2Fl0.i»FlC.0.2F3.0«F7.1,F6.1*I5,F8.0*F6.0»IA»F7.0 

C.F7.1 ,F6.0,F12.6» 

IFC---.6C.0 )"Pfl 

155  CALL  8 ATO  ISTTPP  .A.AWT.T,  I  T  ,C  IS  .CSTP("««)  «VN0W.B!1l 
IFItP8CC.FC.0ICC  TO  160 
tF(CIS.CT.0S-A*U55l,1553 

1551  e--e**ico. 

DC  15*2  LL-I.ICSTP 

1552  R8ATILLl-RBATlCLU*8-l/8* 

PB  |NT  f 6 ,928  18- 

PB|NT*,"  TOTAL  TI-E-  ".IT."  SECS" 

920  FCRPAT C /2C*92F *****  CYCLE  TER-InatEO  OuE  TO  INSUFFICIENT  BATTERY  C 
CAPACITY  TC  -EET  RAfiCE  aeouiPE-ENrS  ****** /35X24HNEXT  ITERATION  Wit 
CL  USEs./AOXlAEBATTERY  -ASS-  .F5.0I 
CO  TC  81 

1553  IF <CE V.LT .OV-AX 1 156 , 160 

156  Pfl AT --PR 4 T- ,  1  . 

R.PRAT«/<PRATA-l.> 

BA-8-*B*».5 

00  157  L-l.ICSTP 

157  B8AT(L 1-R8AT  1 IL »*B«l/8N 
BFLC-RFLO/R 
BARA-R4R-/R 
K-CT-k-CT/R 
A-AX«4-*X*R 
PBlNTT6.925IPBAT-.8- 

PRINT*,"  TOTAL  T I -6 •  -.IT,"  SECS" 

925  FCB *  * T I  25 X 8 1 H « « « »♦  CYCLE  TER-IN.  TEO  OUE  TO  INSUFFICIENT  POWER  TO  - 
CEET  PCvEB  REOLIRE-E.NTS  •  •♦•*, /35X2RHN EXT  ITERATION  WILL  USE5/60X13 
CH-OTCB  PCNEB-  ,  F  *  • l/S0Xl*h8ATTERY  -ASS-  ,F5.0» 

CO  TO  83 
160  CONTINUE 

170  PRIM*,-  TOTAL  CYCLE  TIAE  TPINJ 
TI-E*T«IIT-n/60. 

PR INT I  6 .919)1  1  AE 
919  FCR'ATTFIO.A) 

PRIM*,-  TOTAL  CYCLE  DISTANCE  C K -I »  -• 

PR  INT 16,919  10 
W--V-VM 

PRIM*,"  IOTAL  -OTOR  INPUT  ENERCr  TKW-HR)  -• 

PB  1M  I6.9191WP 

PRINT*,"  TOTAL  Ef'C  IN6  INPUT  ENERCY  IKW-HR)  ■  • 

PBINTI6.919IWP 
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>V*  .vv  /■  .*•  ,%  ,VV  ."O  .>  L>l>  .Y.vjf.  V  ■ 


V.VA  A  A 


PRtNT*,-  TOTAL  REGENERATION  £NE«CY  AVAILASLE  IKW-HRI  .* 
PR  IMIJ.BII  tii« 

3*3  PRINT*. *  8  A  T  T  £  B  Y  CISCHARCE  AT  EM)  OF  CYCLE  It)  •• 

P  8  I  M  <  6  ,  <J  1  9  )  C  !  S 
CAS*WH«.56*f CCST/O 
ELECT. W**EC05T/. 6/0 

PRIM*,*  A V C  CYCLE  ELECTRICAL  COST  IWKH1  -• 

3)0  PR  !NM6,<»l<nCLECT 

PRINT*,*  A VC  FETRCLEUP  FUEL  COST  I4/KNI  •* 

P«  [ST «6,9iq )C AS 
IF  I  IPROC.EO.O  ICO  TO  172 
PRINT* ,*  FINAL  VEHICLE  PARAMETERS  ARE:- 
255  PRINT*, •  PCTCR  POWER*  ".PBatp."  BATTERT  PASS*  ",8* 

PRINT*, •  1CTAL  PASS*  ",P,"  FRONTAL  area*  * »f a 
172  ENO 


1  SUBROUTINE  PTBSPO(YVA,RB,AR,CR,UPSHFT,ONSHFT,I5,K,RPSI 

0  I  .*£  NS  l  ON  GBR  1 c I 

c 

C  this  SUBROUTINE  COMPUTES  POTCR  SPEEO  FOR  PULTI-SPEEO  TRANS* I SS IONS 
5  C  ASC  0»IV|NC  CYCLE  RECUIREPENTS, 

c 


100 

R*S"VVA*6.)*9/fiR*AR*GR(K| 

IF  IRAS.CT.C'PSHFTI  101,103 

101 

IFIK.LT.ISI 102,106 

1C  102 

K  »K  ♦  l 

CC  TO  100 

103 

IF  <R*S.LT.CN3HFTI10*.1C6 

10* 

IFIK.CT.1I  10), 106 

10) 

K-K-l 

1) 

CC  TO  100 

106 

return 

ENO 
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50BBCLMNE  VCt-CP<'<T>P,aAJ.,BfL0,Rp«CT,A»i*,?an«,i<«STP,aPST?,ETA4 
CB8AT,VPAT,IT-*?,N8t=TAV,VHt?,,».S,v,A,VLT.PAvL,PB0CH| 

0  I  -  €  N  S  t  C  N  B *S 1 P  I  10 >  ,E T A« < 10 , 10  ) 

c 

c  1M*  susbcutine  calculates  "crna  po*eb  cevelopeo  anc  cubbent 

C  SfOUfPEfT  FOB  SC«  C.-iCPPfS  vnLTACE  Cu  *•  T  9  CL  L  E  a  S  v[Tm  OC  POTCBS.  IT 

C  A  S  SL " £  5  1.5  VCLTS  JUNCTION  LCSS  ANO  3Z  CCPPEB  LOSS. 

c 

KM«0 

C»B»SABK**OT 

paoc-p»iccc. 

CALL  PTREFFtPBCO.PBATP.B-S.BNSTP.IPSTP.ETAN.ETAl 
IEI PTtP.EC.C I  a  C l , A  02 

A  C I  E  A  I A«PPQO/ ETA*  100. *C/( BFLO«BAflN»BFLO/l BELO-C  1  t 
IP4EAIA.lt.  prcoeau-pbco 
VLT«1£ At  A*#  ARP*RELD*BFLC/C/IRFLD-C1  i  **.5 
VPA».A'P4X*«iaP«RfL0/<RFt0-Cl 
A-E A  lAABFLO/VLT/C 
CO  TO  403 

402  RTOT*BABN»flPLC 

EA!a-PBCC»C/ETA/(C»BTOT) 

IP4EAIA.LT.PPCCIEAIA.PB00 

A«(EA|A/C)*4.S 

V»A*.A*AX«4B8AT«l.03*CaT0T»C>l*1.5 

VLT.A4VP.4X/4PAX 

4C3  IF1VLT.CT.VPAX |KK«1  tVLT«VPA*  1 4  « A  P AX 
IF  4  VLT  .LC . ve  A  1  1404, 406 

404  rFIBK.EO. 01405. 402 

405  PAVL«P 
RETURN 

406  VLT.V8AT 

407  IPfPTTP.EO.C 14C8,*09 

408  EAtA*C*VLT«VLT/RARP/RFLO«4t.-C/RPLOI 
A-EAIA/C/VLIABELO 

CO  TO  410 

405  EAr**lVLT-l.51/(R8AT*l.C3*4RTOT»Cll**2.*C 
A.<£AIA/C1»«.« 

410  CALL  PTBEFP (EA  1  A  ,  P R A T •  ,  B P S . B PS  TP  ,  I  4$ T P , E T A 4 , E T A T 
tP<PTYP. EC. 01411, 412 

411  P A VL  *E  T  A  *|  E  A  I  A ♦ A*A«B  AB VLT* VLT/RFL01 *100000. 

412  PAVL.ETA*(EAIA«A*A«BT0T1/10CCC0. 

CO  TO  413 

413  IPIPAVL.CT.EA  |A/1C00.!PAVL«EA1 A/1C0C. 

IPUHT8.  EC.  1.  ANO.  V.CE.VM  1414,417 

414  PAVLP.f.CCOe7*PP5*HR..<l<n«ETAH/100. 

PBCCH-P-PAVL 

IPIPAVLh.CE.PBCCH 1415, 416 

415  PAVL-P 
RETURN 

416  pavl*pavlh*pavl 

PRQCV.PAVLH 

417  RETURN 


S  U  8  R  C  L  T  f  N  E  YST£PtP,R-S.RAa-,a.<-aT.vBAT,tCSTP,CSTP,A,PAVL.APAX,PRAT 
C^t8»^TPtjn^tP  ,ETA.h,n.R94T  ,  v  L  T  ilhHif*  .ETAH.PRCCH.V.VH.NTYP.RFLOl 
0I*fNSICN  CSTP<5>»RhST?11C1.ETA-110.)01.«BATI5» 

THU  SUPRCUTIFE  CALCULATES  l-CTCR  PCwER  CEVELOPEO  A  N  C  CURRENT 
8ECUIRE0  TOP  FINITE  STEP  VOLTACE  CCNTRCLLERS  W I TH  OC  PCT0«S. 

*••0 

IF  I PTVP.EO.C. AnO.8NS.LT. ICCC. >«HS» 1000. 

501  P*N«1 
C«OPS*RRP0T 
0V*a»R8AT<PI 

5011  VLT«CSTP1N)«VE IT-OV 
1F1PTYP.E0. 05111*5112 

5111  EAIA«C4VLT*VlT/RABN/RFLC«11.-C/8FL0I 
4*EHA*RFL0/VLT/C 

CO  TO  5113 

5112  810T*RFL0*8ARP 

£4  I  *•  t  ivlT-YU  *C/ 1 8  TOT »cn/H  TOT  I  **2.«C 
4«(E4(4/C)*«.5 

5113  OV2«A«B8AT<NI 
lFf4PS(0V2-CV».CT.O.3J5012.5013 

5012  0V-0Y2 

CO  TO  SC  1 1 

5013  CALL  -TPEFF IE  A  1  A ,PB AT*.R*S.R«STP, I 1STP .ETAN.ETAI 
1F1  HYP. EC. 05211 .5212 

5211  P4VL-ET4*(Eil4*4«4«8  4R*..VLT4VLT/8FL0l/100COO. 

CO  TC  5213 

5212  PAVL*£TA«IEAIA»A**2.»flTOT)/lCCCOO. 

5213  IETP4VL.CT.EA14/1C00.IP4VL *£414/1000. 

5015  IF  I •c.FC.Ol 5C2  .507 

502  !F|P,CT.l.4f‘0.4.CT.4"AVI5Cfc«505 

504  1 F I  P4VL.CE . IP  4 .8 1  ) 506 l .505 
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